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Abstract 

The cusp anomalous dimension is a ubiquitous quantity in four-dimensional gauge theories, 
ranging from QCD to maximally supersymmetric A/" = 4 Yang-Mills theory, and it is one of the 
best investigated observables in the AdS/CFT correspondence. In planar A" = 4 SYM theory, 
its perturbative expansion at weak coupling has a finite radius of convergence while at strong 
coupling it admits an expansion in inverse powers of the 't Hooft coupling which is given by a 
non-Borel summable asymptotic series. We study the cusp anomalous dimension in the transition 
regime from strong to weak coupling and argue that the transition is driven by nonperturbative, 
exponentially suppressed corrections. To compute these corrections, we revisit the calculation of 
the cusp anomalous dimension in planar A/" = 4 SYM theory and extend the previous analysis 
by taking into account nonperturbative effects. We demonstrate that the scale parameterizing 
nonperturbative corrections coincides with the mass gap of the two-dimensional bosonic 0(6) 
sigma model embedded into the AdSs x S^ string theory. This result is in agreement with the 
prediction coming from the string theory consideration. 
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1 Introduction 



The AdS/CFT correspondence provides a powerful framework for studying maximally super- 
symmetric M = 4: Yang-Mills theory (SYM) at strong coupling [Ij. At present, one of the best 
studied examples of the conjectured gauge/string duality is the relationship between anomalous 
dimensions of Wilson operators in planar A/" = 4 theory in the so-called SL{2) sector and energy 
spectrum of folded strings spinning on AdSs x [21 [3]. The Wilson operators in this sector are 
given by single trace operators built from L copies of the same complex scalar field and light- 
cone components of the covariant derivatives. These quantum numbers define, correspondingly, 
the twist and the Lorentz spin of the Wilson operators in A/" = 4 SYM theory (for a review, see 
[1]). In dual string theory description [2], [3] they are identified as angular momenta of the string 
spinning on and AdSs part of the background. 

In general, anomalous dimensions in planar A/" = 4 theory in the SL{2) sector are nontrivial 
functions of 't Hooft coupling g'^ = gy^Nc/ {4:7iy and quantum numbers of Wilson operators - 
twist L and Lorentz spin A^. Significant simplification occurs in the limit [5] when the Lorentz 
spin grows exponentially with the twist, L ~ InA^ with A^ —>■ oo. In this limit, the anoma- 
lous dimensions scale logarithmically with A^ for arbitrary coupling and the minimal anomalous 
dimension has the following scaling behavior O El [3, El E] 

77v,l(^/) = [2reusp(^?) + eig,j)] In AT + ... , (LI) 

where j = L/ InN is an appropriate scaling variable and ellipses denote terms suppressed by 
powers of 1/L. Here, the coefficient in front of InA^ is split into the sum of two functions in 
such a way that e{g,j) carries the dependence on the twist and it vanishes for j = 0. The first 
term inside the square brackets in ( 11. ip has a universal, twist independent form [TOl [TT] . It 
involves the function of the coupling constant known as the cusp anomalous dimension. This 
anomalous dimension was introduced in \10\ to describe specific (cusp) ultraviolet divergences of 
Wilson loops |T2| [13] with a light-like cusp on the integration contour p^. The cusp anomalous 
dimension plays a distinguished role in A/" = 4 theory and, in general, in four- dimensional Yang- 
Mills theories since, aside from logarithmic scaling of the anomalous dimension (II. ip . it also 
controls infrared divergences of scattering amplitudes [15], Sudakov asymptotics of elastic form 
factors [16], gluon Regge trajectories [E] etc. 

According to (11. ip . asymptotic behavior of the minimal anomalous dimension is determined 
by two independent functions, rcusp(5') and e{g,j). At weak coupling, these functions are given 
by series in powers of g^ and the first few terms of the expansion can be computed in perturba- 
tion theory. At strong coupling, the AdS/CFT correspondence allows us to obtain expansion of 
rcusp(5') and e{g,j) in powers of 1/g from the semiclassical expansion of the energy of the folded 
spinning string. Being combined together, the weak and strong coupling expansions define asymp- 
totic behavior of these functions at the boundaries of (semi-infinite) interval < g < oo. The 
following questions arise: What are the corresponding interpolating functions for arbitrary g7 
How does the transition from the weak to strong coupling regimes occur? These are the questions 
that we address in this paper. 

At weak coupling, the functions rcusp(fi') and e{g, j) can be found in a generic (supersymmetric) 
Yang-Mills theory in the planar limit by making use of the remarkable property of integrability. 
The Bethe Ansatz approach to computing these functions at weak coupling was developed in 
[T8l [5l [11]. It was extended in [71 [19] to all loops in A/" = 4 SYM theory leading to integral 
BES/FRS equations for rcusp(5') and £{g,j) valid in the planar limit for arbitrary values of the 
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scaling parameter j and the coupling constant g. For the cusp anomalous dimension, the solution 
to the BES equation at weak coupling is in agreement with the most advanced explicit four-loop 
perturbative calculation [20] and it yields a perturbative series for rcusp(5') which has a finite 
radius of convergence [I9] . The BES equation was also analyzed at strong coupling [211 ISS, [23|, [2l] 
but constructing its solution for T cuspid) turned out to be a nontrivial task. 

The problem was solved in Refs. [25l EE] , where the cusp anomalous dimension was found in 
the form of an asymptotic series in 1/g. It turned out that the coefficients of this expansion have 
the same sign and grow factorially at higher orders. As a result, the asymptotic 1/g expansion 
of rcusp(5') is given by a non-Borel summable series which suffers from ambiguities that are 
exponentially small for g ^ oo. This suggests that the cusp anomalous dimension receives 
nonperturbative corrections at strong coupling [25] 

CO 

rcusp(c/) = Yl ^^/^^ ~ ^ "^^'^P + '^("^cusp) ■ (1-2) 

Here the dependence of the nonperturbative scale m^^^^ on the coupling constant rricusp ~ 
gi/i^-ng fQiiQ-yyg^ through a standard analysis [23 EB], from the large order behavior of the 
expansion coefficients, ~ T{k + |) for k ^ oo. The value of the coefficient a in (11. 2p depends 
on the regularization of Borel singularities in the perturbative 1/g expansion and the numerical 
prefactor was introduced for the later convenience. 

Notice that the expression for the nonperturbative scale m^„gp looks similar to that for the 
mass gap in an asymptotically free field theory with the coupling constant ^ 1/g. An important 
difference is, however, that m^^^p is a dimensionless function of the 't Hooft coupling. This is 
perfectly consistent with the fact that Af = 4 model is a conformal field theory and, therefore, 
it does not involve any dimensionfull scale. Nevertheless, as we will show in this paper, the 
nonperturbative scale ml^^^ is indeed related to the mass gap in the two-dimensional bosonic 
0(6) sigma- model. 

The relation (11.21) sheds light on the properties of Tcusplfi') in the transition region (yf ~ 1. 
Going from ^ 1 to (7 = 1, we find that m^^gp increases and, as a consequence, nonperturbative 
O(m^^gp) corrections to Fcuspls') become comparable with perturbative 0(1/5') corrections. We 
will argue in this paper that the nonperturbative corrections play a crucial role in the transition 
from the strong to weak coupling regime. To describe the transition, we present a simplified model 
for the cusp anomalous dimension. This model correctly captures the properties of Tcuspio) at 
strong coupling and, most importantly, it allows us to obtain a closed expression for the cusp 
anomalous dimension which turns out to be remarkably close to the exact value of Fcuspls') 
throughout the entire range of the coupling constant. 

In the AdS/CFT correspondence, the relation (ll.2p should follow from the semiclassical ex- 
pansion of the energy of quantized folded spinning string [21 [3] • In the right-hand side of (II. 2p , 
the coefficient c_i corresponds to the classical energy and describes {k + l)th loop correc- 
tion. Indeed, the explicit two- loop stringy calculation [29] yields the expressions for c_i, cq and 
Ci which are in a perfect agreement with (ll.2p FI However, the semiclassical approach does not 
allow us to calculate nonperturbative corrections to T^nspig) and verification of (II. 2p remains a 
challenge for the string theory. 

Recently, Alday and Maldacena [H] put forward an interesting proposal that the scaling func- 
tion e{g,j) entering (11.10 can be found exactly at strong coupling in terms of nonlinear 0(6) 



^The same result was obtained using different approacli from the quantum string Bcthc Ansatz in Refs. [9[|30|. 
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bosonic sigma model embedded into AdSs x model. More precisely, using the dual description 
of Wilson operators as folded strings spinning on AdSs x and taking into account the one-loop 
stringy corrections to these states [8J, they conjectured that the scaling function e{g,j) should 
be related at strong coupling to the energy density eo{6) in the ground state of the 0(6) model 
corresponding to the particle density po(6) = j/2 



eo(6) 



mo(6) = V/'e--^ [l + 0{l/g)] . 



[1.3) 



This relation should hold at strong coupling and j/mo(6) = fixed. Here the scale mo(6) is 
identified as the dynamically generated mass gap in the 0(6) model with k = 2^/^7i^/^/T{l) 
being the normalization factor. 

The 0(6) sigma model is an exactly solvable theory [31], |32l [331 El] and the dependence of 
eo{6) on the mass scale mo(6) and the density of particles po(6) can be found exactly with a help 
of thermodynamical Bethe ansatz equations. Together with f 1 1.3 1) this allows us to determine 
the scaling function e{g,j) at strong coupling. In particular, for j7mo(6) 1, the asymptotic 
behavior of e{g,j) follows from the known expression for the energy density of the 0(6) model 
in the (nonperturbative) regime of small density of particles [311 El ESI ES] 



9)+j= 



TC 



m 



J 



24 V m 



+ O (f/m^) 



:i.4) 



with m = mo(6)- For j/mo(6) ^ 1, the scaling function e{g,j) admits a perturbative expansion 
in inverse powers of g with the coefficients enhanced by powers of Ini (with i = j / {Ag) ^ 1) [8l[6] 



e{g,j)+j=2f 



1 /3 , 

TT \ 4 



Aix'^g 



*^^-31n. 



+ 4(ln£)2)+0(l//) 



+ 0{t). (1.5) 



This expansion was derived both in string theory [37] and in gauge theory [301 EH E3 yielding 
however different results for the constant go2- The reason for the disagreement remains unclear. 

Remarkably enough, the relation (11.31) was established in planar Af = 4 SYM theory at 
strong coupling [35] using the conjectured integrability of the dilatation operator [7]. The mass 
scale mo(6) was computed both numerically [36J and analytically [351 EH] and it was found to 
be in a perfect agreement with (11. 3p . This result is an extremely nontrivial given the fact that 
the scale mo(6) has a different origin in gauge and in string theory sides of the AdS/CFT. In 
string theory, it is generated by the dimensional transmutation mechanism in two-dimensional 
effective theory describing dynamics of massless modes in the AdSs x sigma model. In gauge 
theory, the same scale parameterizes nonperturbative corrections to anomalous dimensions in 
four- dimensional Yang-Mills theory at strong coupling. It is interesting to note that similar 
phenomenon, when two different quantities computed in four- dimensional gauge theory and in 
dual two-dimensional sigma model coincide, has already been observed in the BPS spectrum in 
M = 2 supersymmetric Yang- Mills theory [lOl HI] . We would like to mention that the precise 
matching of the leading coefficients in perturbative expansion of spinning string energy and 
anomalous dimensions on gauge side was previously found in Refs. H3l B4] . The relation 
( 11.31) implies that for the anomalous dimensions (11.11) the gauge/string correspondence holds at 
the level of nonperturbative corrections. 

As we just explained, the functions rcusp(5') and e{g,j) entering (11.11) receive nonperturbative 
contributions at strong coupling described by the scales mcusp and mo(6), respectively. In A/" = 4 
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SYM theory, these functions satisfy two different integral equations [3, [19] and there is no a 
priori reason why the scales mcusp and r?7,o(6) should be related to each other. Nevertheless, 
examining their leading order expressions, Eqs. (11.21) and fll.3l) . we notice that they have the same 
dependence on the coupling constant. One may wonder whether subleading 0{1/ g) corrections 
are also related to each other. In this paper, we show that the two scales coincide at strong 
coupling to any order oil/ g expansion 

"^cusp = "^0(6) , (1-6) 

thus proving that nonperturbative corrections to the cusp anomalous dimension ( 11.20 and to the 
scaling function (11. 4p are parameterized by the same scale. 

The relations (11.21) and (II. 6p also have an interpretation in string theory. The cusp anomalous 
dimension has the meaning of the energy density of a folded string spinning on AdSs [2|, [6]. As 
such, it receives quantum corrections from both massive and massless excitations of this string 
in the AdSs x sigma model. The 0(6) model emerges in this context as the effective theory 
describing the dynamics of massless modes. In distinction with the scaling function e{g,j), for 
which the massive modes decouple in the limit j/mo(6) = fixed and g ^ oo, the cusp anomalous 
dimension is not described entirely by the 0(6) model. Nevertheless, it is expected that the 
leading nonperturbative corrections to Tcuspig) should originate from nontrivial infrared dynamics 
of the massless excitations and, therefore, they should be related to nonperturbative corrections 
to the vacuum energy density in the 0(6) model. As a consequence, rcusp(5') should receive 
exponentially suppressed corrections proportional to square of the 0(6) mass gap ~ ''^0(6) • 
show in this paper by explicit calculation that this is indeed the case. 

The paper is organized as follows. In Section 2, we revisit the calculation of the cusp anoma- 
lous dimension in planar A/" = 4 SYM theory and construct the exact solution for Fcusplf?)- In 
Section 3, we analyze the obtained expressions at strong coupling and identify nonperturbative 
corrections to Tcnspig)- In Section 4, we compute subleading corrections to the nonperturbative 
scales mcusp and mo(6) and show that they are the same for the two scales. Then, we extend our 
analysis to higher orders in 1/g and demonstrate that the two scales coincide. Section 5 contains 
concluding remarks. Some technical details of our calculations are presented in Appendices. 

2 Cusp anomalous dimension in AT = 4 SYM 

The cusp anomalous dimension can be found in planar Af = 4 SYM theory for arbitrary coupling 
as solution to the BES equation [19j. At strong coupling, Fcusplfi') was constructed in [251 ES] 
in the form of perturbative expansion in 1/g. The coefficients of this series grow factorially 
at higher orders thus indicating that Fcuspls') receives nonperturbative corrections which are 
exponentially small at strong coupling, Eq. (II. 2p . To identity such corrections, we revisit in this 
section the calculation of the cusp anomalous dimension and construct the exact solution to the 
BES equation for arbitrary couphng. 

2.1 Integral equation and mass scale 

In the Bethe ansatz approach, the cusp anomalous dimension is determined by the behavior 
around the origin of the auxiliary function •y^t) related to density of Bethe roots 

rcusp(^?) = -82(?'lim7(t)/t. (2.1) 
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The function 7(t) depends on 't Hooft coupling and has the form 

7(t) = 7+(t)+z7_(t), 



(2.2) 



where 7-|-(t) are real functions of t with a definite parity 7^(±t) = ±7_|_(t). For arbitrary coupling, 
the functions 7-|-(t) satisfy the (infinite-dimensional) system of integral equations 



°° dt 
T 
~ dt 

T 



J2n{t) 



7-(t) 



+ 



7+(t) 

1 _ e-*/(29) ' e*/(25) -1 
7+(t) 7-(t) 



(2.3) 



l_e-*/(29) e*/(29)_i 



with n > 1 and J„(t) being the Bessel functions. These relations are equivalent to BES equa- 
tion [191 provided that 7±(t) verify certain analyticity conditions specified below in Sect. 2.2. 

As was shown in [251 ES] , the equations (12.31) can be significantly simplified with a help of the 
transformation 7(t) r(t): 



r{t) = ( i + zcoth^ )7(t) = r+(t) + ir_(t). 



(2.4) 



We find from (12.11) and (12.41) the following representation for the cusp anomalous dimension 

rcusp(^7) = ~2gT{0) . (2.5) 

It follows from (12. 2p and (12.31) that T^{t) are real functions with a definite parity, r_|_(— t) = 
±Tj_(t), satisfying the system of integral equations 



dt cos(-ut) 
dt sin{ut) 



T_it)-T4t) 
T_it)+T4t) 



: 2 

0, 



(2.6) 



with u being arbitrary real parameter such that — 1 < m < 1. Since T^{t) take real values, we 
can rewrite these relations in a compact form 



dt 



e*"*r_(t) -e-™*r+(t) 



2. 



(2.7) 



To recover (12. 3p . we apply (12. 4p . replace in (12.60 trigonometric functions by their Bessel series 
expansions 



cos u 



sm u 



t) = 2 J^(2n - 1 

n>l 
t) =25^(2 



^ COs((2n - l)(p) J2n-l{t) 
COS if t 

sin(2rav?) J2n{t) 



n>l 



COS 1^ 



t 



{2.i 



^With a slight abuse of notations, we use here the same notation as for Euler gamma-function. 
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with u = sinyj, and finally compare coefficients in front of cos ((2n — l)ip)/cos(p and sin{2nip)/cosip 
in both sides of (12. 6p . It is important to stress that, doing this calculation, we interchanged the 
sum over n with the integral over t. This is only justified for (p real and, therefore, the relation 
(EH) only holds for -1 < m < 1. 

Comparing (12. 7p and (12.31) we observe that the transformation 7^ eliminates the 

dependence of the integral kernel in the left-hand side of ( 12. 7p on the coupling constant. One 
may then wonder where does the dependence of the functions r_|_(t) on the coupling constant 
come from? We will show in the next subsection that it is dictated by additional conditions 
imposed on analytical properties of solutions to (12.70 . 

The relations (12.50 and (12.60 were used in [25j to derive asymptotic (perturbative) expansion of 
rcusp(5) in powers of 1/g. This series suffers however from Borel singularities and we expect that 
the cusp anomalous dimension should receive nonperturbative corrections ~ e~'^'^3 exponentially 
small at strong coupling. As was already mentioned in the Introduction, similar corrections are 
also present in the scaling function e{g,j) which controls asymptotic behavior of the anomalous 
dimensions (II. ip in the limit when Lorentz spin of Wilson operators grows exponentially with 
their twist. According to (11.30 . for j7mo(6) = fixed and g ^ 00, the scaling function coincides 
with the energy density of the 0(6) model embedded into AdSs x S^. The mass gap of this model 
defines a new nonperturbative scale mo(6) in the AdS/CFT. Its dependence on the coupling g 
follows univocally from the FRS equation and it has the following form [351 EH] 

mo(6) = ^e--^-^e--^Re / [V ^{t) + _{t)) , (2.9) 

where T_^{t) are solutions to (12. 7p . To compute the mass gap (12. 9p . we have to solve the integral 
equation (12. 7p and, then, substitute the resulting expression for T_^(t) into (12.90 . Notice that the 
same functions also determine the cusp anomalous dimension (12. 5p . 

Later in the paper, we will construct the solution to the integral equation (12.70 and, then, 
apply (12. 5p to compute nonperturbative corrections to Fcusp(fi') at strong coupling. 



2.2 Analyticity conditions 

The integral equations (12.70 and (12.30 determine T_^{t) and 'y±it), or equivalently the functions 
F(t) and 7(t), up to a contribution of zero modes. The latter satisfy the same integral equations 
(12.71) and (12. 3p but without inhomogeneous term in the right-hand side. 

To fix the zero modes, we have to impose additional conditions on solutions to (12. 7p and (12. 3p . 
These conditions follow unambiguously from the BES equation [23l [25] and they can be formu- 
lated as a requirement that 7-|-(t) should be entire functions of t which admit a representation in 
the form of Neumann series over Bessel functions 

7_(t) = 2 ^ (2n - 1) J2n-l(t)72n-l , (2.10) 
7+(t) =2^ {2n) J2n{t)l2n, 

with the expansion coefficients 72^-1 and 72^ depending on the coupling constant. This implies 
in particular that the series on the right-hand side of (I2.10p are convergent on the real axis. Using 
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orthogonality conditions for the Bessel functions, we obtain from (]2.10p 

/"°° dt dt 

72„-l= / -J2n-l(t)7-W, 72n= / -^2n(t)7+W- (2.11) 

Jo Jo 

Here we assumed that the sum over n in the right-hand side of (12.101) can be interchanged with 
the integral over t. We will show below that the relations fl2.10p and (12.111) determine a unique 
solution to the system (12.31) . 

The coefficient 71 plays a special role in our analysis since it determines the cusp anomalous 
dimension (12. ip . 

rcusp(^?) = 8/71(^7)- (2.12) 

Here we applied (12.20 and (12.1 OP and took into account small— t behavior of the Bessel functions, 
Jnit) ust^O. 

Let us now translate (12.100 and (12.110 into properties of the functions r^(t), or equivalently 
r(t). It is convenient to rewrite the relation (12.40 as 

= ^^^^^ sin(f)sin(f) = n t-4.,fc " ^''''^ 

^^9' k=-oo ^ 

Since 7 (it) is an entire function in the complex t— plane, we conclude from (12.130 that T{it) has 
an infinite number of zeros, T{it^cros) = 0, and poles, T{it) ~ l/(t — tpoies), on real t— axis located 
at 

tzeros = 47r^ " i) , tpoles = 47r^f , (2.14) 

where £, G Z and ^' 7^ so that T{it) is regular at the origin (see Eq. (12. ip ). Notice that T(it) 
has an additional (infinite) set of zeros coming from the function 7 (it) but, in distinction with 
(I2.14p . their position is not fixed. Later in the paper we will construct the solution to the integral 
equation (12. 6p which satisfies the relations (I2.14p . 

2.3 Toy model 

To understand the relationship between analytical properties of T{it) and properties of the cusp 
anomalous dimension, it is instructive to slightly simplify the problem and consider a 'toy' model 
in which the function r(zt) is replaced with V'^^°^\it). 

We require that r(*°y)(it) satisfies the same integral equation (12. 6p and define, following (12. 5p . 
the cusp anomalous dimension in the toy model as 

rS(^?) = -2^7r(*°yno). (2.15) 

The only difference compared to V{it) is that V^'^°^\it) has different analytical properties dictated 
by the relation 

r(t°y)(^t) =7(*°y)(zt)^^, (2.16) 



while 7'^*°y^(zt) has the same analytical properties as the function 7(it)|j This relation can be 
considered as a simplified version of (12.130 . Indeed, it can be obtained from (12.130 if we retained 



•^Notice that the function ^'^'^°^\t) does not satisfy the integral equation (|2.3p anymore. Substitution of (|2.16p 
into (|2.7p yields integral equation for 7*^*°y^(t) which can be obtained from (|2.3p by replacing l/(l — e^*/^^^)) 
^ + i and l/(e*/(29) _i) ^ ££ _ 1 in the kernel in the left-hand side of ([Qjl . 



8 



in the product only one term with k = 0. As compared with (I2.14p . the function T^^'^^^it) does 
not have poles and it vanishes for t = —ng. 

The main advantage of the toy model is that, as we will show in Sect. 12. 8[ the expression for 
Tcu'sp ((?) can be found in a closed form for arbitrary value of the coupling constant (see Eq. (12.501) 
below). We will then compare it with the exact expression for T cuspid) and identify the difference 
between the two functions. 



2.4 Exact bounds and unicity of the solution 

Before we turn to finding the solution to (12.61) . let us demonstrate that this integral equation 
supplemented with the additional conditions ( I2.10p and ( 12. lip on its solutions, leads to nontrivial 
constraints for the cusp anomalous dimension valid for arbitrary coupling g. 

Let us multiply both sides of the two relations in (12.31) by 2{2n — l)72n-i and 2(2?7,)72„, 
respectively, and perform summation over n > 1. Then, we convert the sums into the functions 
7_i_(t) using (12.101) and add the second relation to the first one to obtain 

^' = L T l-e-V(^.) • ^^-^'^ 

Since 7_|_(t) are real functions of t and the denominator is positively definite for < t < oo, this 
relation leads to the following inequality 



7i 



> 



jyj{l.{t)f>2^l>Q. (2.18) 



Here we replaced the function 7_. {t) by its Bessel series (I2.10p and made use of the orthogonality 
condition for the Bessel functions with odd indices. We deduce from (I2.18P that 

< 71 < i (2.19) 

and, then, apply (12.120 to translate this inequality into the following relation for the cusp anoma- 
lous dimension 

0<reusp((7) < V. (2.20) 

We would like to stress that this relation should hold in planar Af = 4 SYM theory for arbitrary 
coupling g. 

Notice that the lower bound on the cusp anomalous dimension, Fcuspls') > 0, holds in any 
gauge theory [llj. It is the upper bound Tcuspig) < 4(yf^ that is a distinguished feature of A/" = 4 
theory. Let us verify the validity of (I2.20p . At weak coupling Fcuspls') admits perturbative 
expansion in powers of g^ [20] 



rcusp(5') = % 



3 " 45 V630 
while at strong coupling it has the form [251 EHl [26] 



rcusp(5') = 25- 



31n2 _i K _2 /3Kln2 27(3 \ s ^, _4^ 



-IM^'' -1-64^ + 2048;^'^" +^^^ 



(2.21) 



(2.22) 



*Our analysis here goes along the same lines as in Appendix A of [35]. 
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with K being the Catalan constant. It is easy to see that the relations (]2.2ip and (12.22^ are in 
an agreement with (12.201) . 

For arbitrary g we can verify the relation (I2.20p by using the results for the cusp anomalous 
dimension obtained from numerical solution of the BES equation [251 US]- The comparison is 
shown in Figure 1. We observe that the upper bound condition Tcnspid)/ C^g) < '^9 is indeed 
satisfied for arbitrary g > 0. 




Fi gure 1: Dependence of the cusp anomalous dimension Tcuspig)/ i'^g) on the coupling constant. Dashed 
Hne denotes the upper bound 2g. 

We are ready to show that the analyticity conditions formulated in Sect. 2.2 specify a unique 
solution to (12. 3p . As was already mentioned, solutions to (12.31) are defined modulo contribution 
of zero modes, 7(t) — > 7(t) +7^°^(t), with 'y^^\t) being solution to homogenous equations. Going 
through the same steps that led us to (12.171) we obtain 







t 1 - e-*/(25) 



(0), 



(2.23) 



where zero on the left-hand side is due to absence of the inhomogeneous term. Since the integrand 
is a positively definite function, we immediately deduce that 7*-°''(t) = and, therefore, the 
solution for 7(t) is unique. 



2.5 Riemann-Hilbert problem 

Let us now construct the exact solution to the integral equations (12.71) and (12.31) . To this end, it 
is convenient to Fourier transform the functions (12. 2p and (12.40 



f e^^*7(t). 



(2.24) 
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According to ( I2.2p and (I2.10p . the function 7(t) is given by the Neumann series over Bessel 
functions. Then, we perform the Fourier transform on both sides of (12.101) and use the well- 
known fact that the Fourier transform of the Bessel function Jn{t) vanishes for A;^ > 1 to deduce 
that the same is true for 7(t) leading to 

7(fc) = 0, forfc2>l. (2.25) 

This implies that the Fourier integral for 7(t) only involves modes with — 1 < A; < 1 and, 
therefore, the function 7(t) behaves at large (complex) t as 

7(t) ~ el*l , for \t\ oo. (2.26) 

Let us now examine the function T{k). We find from (I2.24p and (I2.13p that T{k) admits the 
following representation 

df , sinh(:r + ^t) 

m= ^e^'' ■ . \ ■ M t). (2.27) 
27r smh{j^) sm(^) 

Here the integrand has poles along the imaginary axis at t = Anign (with n = ±1, ±2, . . .). 

It is suggestive to evaluate the integral (12.271) by deforming the integration contour to infinity 
and by picking up residues at the poles. However, taking into account the relation (12.261) . we find 
that the contribution to (I2.27P at infinity can be neglected for k"^ > 1 only. In this case, closing 
the integration contour into the upper (or lower) half-plane for k > 1 (or k < —1) we find 



r(A;) 9{k ^) e~^™^(^-') + 9{-k - 1) J] 9) e-^^'^'si-k-i) _ (2.28) 

n>l n>l 



Here the notation was introduced for A;— independent expansion coefficients 

c±in, g) = T^g-i{±ATxign) e"^'^"^ , (2.29) 

where the factor e~^™^ is inserted to compensate exponential growth of 7(±47ri(7n) ~ e^™^ at 
large n (see Eq. (I2.26P ). For /c^ < 1, we are not allowed to neglect the contribution to (12.271) at 
infinity and the relation (12.28^ does not hold anymore. As we will see in a moment, for k"^ < 1 
the function T{k) can be found from (12. 7p . 

Comparing the relations (I2.25P and (I2.28p . we conclude that, in distinction with 7(A;), the 
function T{k) does not vanish for k'^ > 1. Moreover, each term in the right-hand side of (12.280 is 
exponentially small at strong coupling and the function scales at large k as T{k) ~ e~^'^^(l^l~^). 
This implies that nonzero values of T{k) for k"^ > 1 are of nonperturbative origin. Indeed, in 
perturbative approach of [2S], the function r(t) is given by the Bessel function series analogous to 
(I2.10p and, similar to (12.251) . the function T{k) vanishes for fc^ > 1 to any order inl/g expansion. 

We note that the sum in the right-hand side of (I2.28P runs over poles of the function T{it) 
specified in (^J^. We recall that in the toy model (I2A6D . T'-^°^\it) and j'-^°y\it) are entire 
functions of t. At large t they have the same asymptotic behavior as the Bessel functions, 
'jy^"y) (it) ~ e^**. Performing their Fourier transformation (I2.24p . we find 

^(toy) _ p(toy) (^) = , for k^>l, (2.30) 



'We recall that j{t) = 0{t) and, therefore, the integrand is regular at t = 0. 
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in a close analogy with (I2.25p . Comparison with (12.28P shows that the coefficients (12.29P vanish 
in the toy model for arbitrary n and g 



.(toy) , 



{n,g) = ct°^\n,g) = 0. 



(2.31) 



The relation (12.281) defines the function T{k) for k"^ > 1 but it involves the coefficients c±(n, g) 
that need to be determined. In addition, we have to construct the same function for fc^ < 1. 
To achieve both goals, let us return to the integral equations (12. 6 p and replace r_|_(t) by Fourier 
integrals (see Eqs. ^TIM and (CT ) 



r+(t) 

T_{t) 



dk cos{kt) T{k) , 

oo 

dk sin(H) r(A;) . 



(2.32) 



In this way, we obtain from (12. 6p the following remarkably simple integral equation for T[k) 

dkf{k) 



k — u 



+ TxTiu) = -2. 



~l<u< 1 



(2.33) 



where the integral is defined using the principal value prescription. This relation is equivalent to 
the functional equation obtained in [26] (see Eq. (55) there). 

Let us split the integral in (I2.33P into k"^ < 1 and k"^ > 1 and rewrite (I2.33P in the form of 
singular integral equation for the function r(A;) on the interval —l<k<l 



, 1 dkT(k) 



71 J_i k — u 
where the inhomogeneous term is given by 

0(u) = --|2 + 

TT 



dkVjk) ^ 

^ k — u 



-1 < M < 1) 



dkVjk) 
k — u 



(2.34) 



(2.35) 



Since integration in (12.351) goes over fc^ > 1, the function T{k) can be replaced in the right-hand 
side of (12.351) by its expression (I2.28P in terms of the coefficients c±{n,g). 

The integral equation (12.340 can be solved by standard methods |16]. A general solution for 
f{k) reads (for -1 < A; < 1) 



v{k) = -m 



l_ (l + k\ '^V' _ V2 ( l + k 



2Tx\l-k 



I u — k \1 + u 



n \1~ k J l + k 



1/4 



(2.36) 



where the last term describes the zero mode contribution with c being an arbitrary function of 
the coupling. We replace (f){u) by its expression (I2.35p . interchange the order of integration and 
find after some algebra 



m 



TT \l-k 



1/4 



C 1 

l + k ^ 2 



=° dpTjp) f p-1 
oo P-k \p+l 



1/4 



(2.37) 
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Notice that the integral in the right-hand side of (I2.37P goes along the real axis except the interval 
[— 1, 1] and, therefore, T{p) can be replaced by its expression (12.281) . 

Being combined together, the relations (12.281) and (I2.37P define the function r(A;) for — oo < 
/c < cxD in terms of (an infinite) set of yet unknown coefficients c±{n,g) and c{g). To fix these 
coefficients we will first perform Fourier transform of r(A;) to obtain the function r{t) and, then, 
require that r(t) should have correct analytical properties (12.141) . 



2.6 General solution 



We are now ready to write down a general expression for the function T{t). According to ( ]2.24p . 
it is related to T{k) through the inverse Fourier transformation 

"1 



-1 



dke 



-ikt ■ 



T{k)+ / dke-'^'V^k) 



(2.38) 



where we split the integral into three terms since T{k) has a different form for k < —1, —1 < k < 1 
and k > 1. Then, we use the obtained expressions for r(A;), Eqs. (12.281) and (I2.37p . to find after 
some algebra the following remarkable relation (see Appendix [B] for details) 



T{it) = ut)v,{t) + h{m{t). 

Here the notation was introduced for 



hit) 

m 



n>l 



Anng — t 



n>l 



c+(n,^)-^^ ^ + c_(n,^) 



^T:ng + 1 

Uoi^Trng) 



(2.39) 



(2.40) 



Anng — t 



Aimg + 1 



Also, Vn and (with n = 0, 1) stand for integrals 



VJx) 



TT 



rfu(l + M)^/^-"(l 



u)-^/^e"" 



(2.41) 



-1 



1 1"°° 

U^{x) = - du {u ± lY^'\u T l)i/4-ne-("-i)" , 

2 Ji 



which can be expressed in terms of Whittaker functions of 1st and 2nd kind [17] (see Appendix [D|) . 
We would like to emphasize that the solution (12.391) is exact for arbitrary coupling g > ^ and 
that the only undetermined ingredients in (12.391) are the expansion coefficients c±{n,g) and c{g). 

In the special case of the toy model, Eq. (I2.3ip . the expansion coefficients vanish, c^°^\n, g) = 
0, and the relation (I2.40p takes a simple form 



/f^)(t) 



.(toy) 



(2.42) 



Substituting these expressions into (I2.39P we obtain a general solution to the integral equation 
(12. 7p in the toy model 

r(t°y)(it) = -Voit) - c^'°y\g)Vi{t) . (2.43) 
It involves an arbitrary (7— dependent constant c*^*"^-* which will be determined in Sect. 12.81 
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2.7 Quantization conditions 



The relation f l2.39p defines a general solution to the integral equation (12 .7^ . It still depends on 
the coefficients c±{n,g) and c{g) that need to be determined. We recall that T{it) should have 
poles and zeros specified in (12.141) . 

Let us first examine poles in the right-hand side of ( 12.390 . It follows from ( 12.4ip that Vq (t) and 
Vi{t) are entire functions of t and, therefore, poles can only come from the functions fo{t) and 
flit). Indeed, the sums entering (I2.40p produce an infinite sequence of poles located at t = ±4™ 
(with n > 1) and, as a result, the solution (I2.39P has a correct pole structure (12.140 . Let us now 
require that T{it) should vanish for t = t^em specified in (I2.14p . This leads to an infinite set of 
relations 

r (47rzc/(£- i)) = 0, ieZ. (2.44) 

Replacing T{it) by its expression (I2.39p . we rewrite these relations in equivalent form 

/o (te) Vo {U) + h{U)Vi{U) = , t, = A7ig{^-\). (2.45) 

The relations (I2.44p and (12.450 provide the quantization conditions for the coefficients c{g) and 
c±{n,g) that we will analyze in Sect. [31 

Let us substitute (12.390 into the expression (12. 5p for the cusp anomalous dimension. The 
result involves the functions Vn{t) and fn{t) (with n = 1,2) evaluated at t = 0. It is easy to see 
from flCTj) that ^(O) = 1 and Vi{0) = 2. In addition, we obtain from flCT]) that /o(0) = -1 
for arbitrary coupling leading to 

reusp(^?) = 2(7[l-2/i(0)]. (2.46) 

Replacing /i(0) by its expression (12.401) we find the following relation for the cusp anomalous 
dimension in terms of the coefficients c and c± 

T,^M = 2g\l + 2c{g)-2Y,[c~{n,g)Uo{4nng) + c+{n,g)U+{4nng^ I . (2.47) 

I n>l J 

We would like to stress that the relations (12.460 and (12.471) are exact and hold for arbitrary 
coupling g. This implies that, at weak coupling, it should reproduce the known expansion of 
Tcuspig) in positive integer powers of g^ [20j. Similarly, at strong coupling, it should reproduce 
the known 1/g expansion [251125] and, most importantly, describe nonperturbative, exponentially 
suppressed corrections to Tcuspig)- 



2.8 Cusp anomalous dimension in the toy model 

As before, the situation simplifies for the toy model (12.430 . In this case, we have only one 
quantization condition T'^^"^\~7rig) = which follows from (12.160 . Together with (12.430 it allows 
us to fix the coefficient c^*"^-* (g) as 



.(toy) 



Vo{-ng) 



V.i-irg) • 

Then, we substitute the relations (I2.48P and ( 12.310 into ( 12.470 and obtain 

M-rrg) 



rS(^) = 2^7 [I + 2c^'°y\g)] = 2g 



'Vii-^g)\ 



(2.48) 



(2.49) 
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Replacing Vo(— tt^j) and Vi(— vr^f) by their expressions in terms of Whittaker function of the first 
kind (see Eq. (ID. 21) ). we find the following remarkable relation 



T^:S{g) = 2g 



1 - {27rg) 



■ 1/2 ^1/4,1/2 (27r5f) 
M_i/4,o(27r^) 



(2.50) 



which defines the cusp anomalous dimension in the toy model for arbitrary coupling g > 0. 

Using fl2.50p it is straightforward to compute riu°sp'(5') for arbitrary positive g. By construc- 
tion, riu°p ((yf) should be different from rcusp(5')- Nevertheless, evaluating (12.501) for < (7 < 3, 
we found that the numerical values of FclSp (fif) are very close to the exact values of the cusp 
anomalous dimension shown by the solid line in Figure 1. Also, as we will show in a moment, 
the two functions have similar properties at strong coupling. To compare these functions, it is 
instructive to examine the asymptotic behavior of rcu°sp''((7) at weak and at strong coupling. 



2.8.1 Weak coupling 

At weak coupling, we find from (12.501) 

rS (9) = 1-9' -I -V - ^ rrY + ^ - ^ - V - ^ + 0(/) • (2.51) 

Comparison with (I2.2ip shows that this expansion is quite different from the weak coupling 
expansion of the cusp anomalous dimension. In distinction with rcusp(fi'), the expansion in (12.511) 
runs both in even and odd powers of the coupling. In addition, the coefficient in front of (7" in 
the right-hand side of (I2.5ip has transcendentality {n — 1) while for Tcnsp{9) it equals [n — 2) 
(with n taking even values only). 

Despite of this and similarly to the weak coupling expansion of the cusp anomalous dimen- 
sion [19], the series (I2.5ip has a finite radius of convergence \go\ = 0.796. It is determined by 
the position of the zero of the Whittaker function closest to the origin, M_i/4 o(27r(yfo) = for 
go = —0.297 ± i 0.739. Moreover, numerical analysis indicates that riu°'p(5') has an infinite num- 
ber of poles in the complex (7— plane. The poles are located in the left-half side of the complex 
plane. Keg < 0, symmetrically with respect to the real axis, and they approach progressively the 
imaginary axis as one goes away from the origin. 



2.8.2 Strong coupling 

At strong coupling, we can replace the Whittaker functions in (I2.50p by their asymptotic expan- 
sion for (7 ^ 1. It is convenient however to apply (I2.49P and replace the functions Vol—Trg) and 
Vi{—7Tg) by their expressions given in (1D.14P and (]D.16p . respectively. In particular, we have 
(see Eq. dHU) 



Vo{-7ig) = e 



l/(2a) ^ 



5/4 



r(|) 



F{ll\a + tO)+A'F{-ll\-a) 



a = l/{2ng), (2.52) 



where the parameter A^ is defined as 



A2 = aa-i/2g-i/aU4 



r(| 



5\ ' 



a 



-3i7r/4 



(2.53) 
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Here, b\ — a) is expressed in terms of the confluent liypergeometric function of the second 
kind (see Eqs. flDTT2l) and (E7j) in Appendix D and Eq. fl236D below) [17] 



F(i,f|-a) =a-^/^[/o+(l/(2«))/r(|), (2.54) 

F(-i,f|-a)=a-^/Vo-(l/(2«))/r(|). 

The function F{a, b\ — a) defined in this way is an analytical function of a with a cut along the 
negative semi- axis. 

For positive a = l/{2TTg), the function F(— i, || — a) entering fl2.52p is defined away from 
the cut and its large g expansion is given by Borel summable asymptotic series (for a = and 

F{a, b\-a) = J2 \^ ^ rJpL ' = 1 - «a6 + 0{a') , (2.55) 

fc>o ■ y ) y ) 

with the expansion coefficients growing factorially to higher orders in a. This series can be 
immediately resummed by means of the Borel resummation method. Namely, replacing r(a + k) 
by its integral representation and performing the sum over k we find for Re a > 

F(a,b\-a) = —- dss"~Vl + s)-^e-"/", (2.56) 

r(a) Jo 

in agreement with (12.541) and (I2.4ip . 

The relation (12.551) holds in fact for arbitrary complex a and the functions F{a,b\a ±«0), 
defined for a > above and below the cut, respectively, are given by the same asymptotic 
expansion (I2.55P with a replaced by —a. The important difference is that now the series (12.550 
is not Borel summable anymore. Indeed, if one attempted to resum this series using the Borel 
summation method, one would immediately find a branch point singularity along the integration 
contour at s = 1 

-a rco 

F(a,b\a±iO) = —- ds s^-^l - s T lO)-'' e-'/'' . (2.57) 

r(a) Jo 

The ambiguity related to the choice of the prescription to integrate over the singularity is known 
as Borel ambiguity. In particular, deforming the s— integration contour above or below the 
cut, one obtains two different functions F{a,b\a ± iO). They define analytical continuation of 
the same function F{a,b\ — a) from Rea > to the upper and lower edge of the cut running 
along the negative semi-axis. Its discontinuity across the cut, F(a, b\a + iO) — F{a, b\a — iO) is 
exponentially suppressed at small a > and is proportional to the nonperturbative scale (see 
Eq. (ID.17P ). This property is perfectly consistent with the fact that the function (12.520 is an 
entire function of a. Indeed, it takes the same form if one used a — iO prescription in the first 
term in the right-hand side of (12.520 and replaced a in (12.531) by its complex conjugated value. 

We can now elucidate the reason for decomposing the entire Vq— function in (I2.52p into the 
sum of two F— functions. In spite of the fact that analytical properties of the former function are 
simpler compared to the latter functions, its asymptotic behavior at large g is more complicated. 
Indeed, the F— functions admit asymptotic expansions in the whole complex g— plane and they 
can be unambiguously defined through the Borel resummation once their analytical properties are 
specified (we recall that the function F{a, b\a) has a cut along positive semi-axis). In distinction 
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with this, the entire function Vo(— vr^f) admits different asymptotic behaviors for positive and 
negative values of g in virtue of the Stokes phenomenon. Not only does it restrict the domain 
of validity of each asymptotic expansion, but it also forces us to keep track of both perturbative 
and nonperturbative contributions in the transition region from positive to negative including 
the transition from the strong to weak coupling. 

We are now in position to discuss the strong coupling expansion of the cusp anomalous dimen- 
sion in the toy model, including into our consideration both perturbative and nonperturbative 
contributions. Substituting f l2.52p and similar relation for Vi{—Tig) (see Eq. ( ]D.16P ) into (12.491) 
we find (for a"*" = a + and a = l/{2iTg)) 

Since the parameter is exponentially suppressed at strong coupling, Eq. (I2.53p . and, at the 
same time, the F— functions are all of the same order, it makes sense to expand the right-hand 
side of (12.581) in powers of and, then, study separately each coefficient function. In this way, 
we identify the leading, A^ independent term as perturbative contribution to rcu°p ((?) and the 
O(A^) term as the leading nonperturbative correction. More precisely, expanding the right-hand 
side of (12.581) in powers of A^ we obtain 

Tt^ig)/{2g) = Coia) - aA'C^ia) + ^a'A^da) + ©(A^) . (2.59) 

Here the expansion runs in even powers of A and the coefficient functions Cfc(a) are given by 
algebraic combinations of F— functions 

Co = l-a f ^ , C2 = 2, C4= \ , (2.60) 



where we applied (]D.9P and (1D.12|) to simplify the last two relations. Since the coefficient 



functions are expressed in terms of the functions F{a,b\a~^) and F{a,b\ — a) having the cut 
along the positive and negative semi-axis, respectively, Ckia) are analytical functions of a in the 
upper-half place. 

Let us now examine the strong coupling expansion of the coefficient functions (12.601) . Replac- 
ing F— functions in (12.601) by their asymptotic series representation (I2.55P we get 

Co=l-a--a'--a^- -a^ - —a' + O (a^) , (2.61) 
4 8 64 128 ^ ^ ' ^ ' 

^ , 1 11 2 151 3 13085 4 ^.5. 
C2 = 1 a a-^ a* + O (a^) , 



Not surprisingly, these expressions inherit the properties of the F— functions - the series (12.610 
are asymptotic and non-Borel summable. If one simply substituted the relations (12.611) into the 
right-hand side of (I2.59p . one would then worry about the meaning of nonperturbative O(A^) 
corrections to (I2.59P given the fact that the strong coupling expansion of perturbative contribution 
Co (a) suffers from Borel ambiguity. We recall that appearance of exponentially suppressed 
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61 4 433 
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-a'^ 
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64 128 
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corrections to rcu°p ((?) is ultimately related to the Stokes phenomenon for the function Vo{—TTg), 
Eq. fl2.52p . As was already mentioned, this does not happen for the F— function and, as a 
consequence, its asymptotic expansion, supplemented with the additional analyticity conditions, 
allows us to reconstruct the F— function through the Borel transformation, Eqs. fl2.56p and f l2.57p . 
Since the coefficient functions fl2.60l) are expressed in terms of the F— functions, we may expect 
that the same should be true for the C— functions. Indeed, it follows from the unicity condition of 
asymptotic expansion [27j, that the functions Co(q;), 6*2(0;), C^i^a), . . . are uniquely determined by 
their series representations fl2.6ip as soon as the latter are understood as asymptotic expansions 
for the functions analytical in the upper half plane Ima > 0. This implies that the exact 
expressions for the functions fl2.60p can be unambiguously constructed by means of the Borel 
resummation but the explicit construction remains beyond the scope of the present study. 

Since the expression fl2.58p is exact for arbitrary coupling g we may now address the question 
formulated in the Introduction: how does the transition from the strong to the weak coupling 
regime occur? We recall that, in the toy model, V'^uip (g) / {2g) is given for -C 1 and g ^ 1 hy 
the relations fl2.5ip and (12.590 . respectively. Let us choose some sufficiently small value of the 
coupling constant, say g = 1/4, and compute Tcuip (g) / {2g) using three different representations. 
Firstly, we substitute g = 0.25 into (12.580 and find the exact value as 0.4424(3). Then, we use 
the weak coupling expansion (I2.5ip and obtain a close value 0.4420(2). Finally, we use the strong 
coupling expansion (I2.59P and evaluate the first few terms in the right-hand side of (12.590 for 
g = 0.25 to get 

Eq. (12391) = (0.2902-0. 1434 i) + (0.1517 + 0.1345 i) 

+ (0.0008 + 0.0086 i) - (0.0002 - 0.0003 i) + . . . = 0.4425 + . . . (2.62) 

Here the four expressions inside the round brackets correspond to contributions proportional to 
A°, A^, A"^ and A^, respectively, with A'^{g = 0.25) = 0.3522 x e"^*'^/^ being the nonperturbative 
scale (12331) . 

We observe that each term in (I2.62p takes complex values and their sum is remarkably close 
to the exact value. In addition, the leading O(A^) nonperturbative correction (the second term) 
is comparable with the perturbative correction (the first term). Moreover, the former term starts 
to dominate over the latter one as we go to smaller values of the coupling constant. Thus, 
the transition from the strong to weak coupling regime is driven by nonperturbative corrections 
parameterized by the scale A^. Moreover, the numerical analysis indicates that the expansion of 
Fcusp ((?) in powers of A^ is convergent for Re (7 > 0. 



2.8.3 Prom toy model to the exact solution 

The relation (I2.59P is remarkably similar to the expected strong coupling expansion of the cusp 
anomalous dimension (11.20 with the function Co(a) providing perturbative contribution and 
A^ defining the leading nonperturbative contribution. Let us compare Co (a) with the known 
perturbative expansion (12.220 of rcusp(5')- In terms of the coupling a = l/{27ig), the first few 
terms of this expansion look as 

, , 31n2 K . /3Kln2 27C3\ , 

r^^M/i^g) = 1 - ^« - - [—^ + ^ J a' + ..., (2.63) 

where ellipses denote both higher order corrections in a and nonperturbative corrections in A^. 
Comparing (I2.63P and the first term, Co{a), in the right-hand side of (I2.59p . we observe that 
both expressions approach the same value 1 as a 0. 
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As was already mentioned, the expansion coefficients of the two series have different transcen- 
dentahty - they are rational for the toy model, Eq. f l2.6ip . and have maximal transcendentality 
for the cusp anomalous dimension, Eq. fl2.63p . Notice that the two series would coincide if one 
formally replaced the transcendental numbers in (12.631) by appropriate rational constants. In 
particular, replacing 

31n2 K 1 9C3 1 , , 

^1, — 2.64 

2 ' 2 2' 32 3 ' ' ^ ^ 

one obtains from (12.631) the ffist few terms of perturbative expansion (12.611) of the function Cq 
in the toy model. This rule can be generalized to all loops as follows. Introducing an auxiliary 
parameter r, we define the generating function for the transcendental numbers in (12.641) and 
rewrite (12.641) as 



exp 



31n2 K , QCs 3^ 

T T H r + . . 

2 2 32 



exp 



r2 

y + y + 



(2.65) 



Going to higher loops, we have to add higher order terms in r to both exponents. In the right- 
hand side, these terms are resummed into exp(ln(l + t)) = 1 + r, while in the left-hand side they 
produce the ratio of Euler gamma-functions leading to 



r(i)r(i + ^)r(|-i) 
r(f)r(i-i)r(i + 5) 



;i + r). (2.66) 



Taking logarithms in both sides of this relation and subsequently expanding them in powers of r, 
we obtain the subtitution rules which generalize (12.641) to the complete family of transcendental 
numbers entering into the strong coupling expansion (I2.63p . At this point, the relation (12.660 
can be thought of as an empirical rule, which allows us to map the strong coupling expansion of 
the cusp anomalous dimension ( 12.63^ into that in the toy model, Eq. ( 12.6ip . We will clarify its 
origin in Sect. 4.2. 

In spite of the fact that the numbers entering both sides of (I2.64p have different transcenden- 
tality, we may compare their numerical values. Taking into account that 3 In 2/2 = 1.0397(2), 
K/2 = 0.4579(8) and 9(3/32 = 0.3380(7) we observe that the relation (12:61) defines a meaningful 
approximation to the transcendental numbers. Moreover, examining the coefficients in front of 
r" in both sides of (I2.65P at large n, we find that the accuracy of approximation increases as 
n — » 00. This is in agreement with the observation made in the beginning of Sect. 2.8, that 
the cusp anomalous dimension in the toy model rcu°sp'(5') is close numerically to the exact ex- 
pression Fcusplf?)- In addition, the same property suggests that the coefficients in the strong 
coupling expansion of riu°sp ((?) and Tcusp{g) should have the same large order behavior. It was 
found in [25j that the expansion coefficients in the right-hand side of (I2.63P grow at higher orders 
as rcusp(fl') ~ ^k^i^ + !)«''• It is straightforward to verify using (I2.60p and (12.550 that the 
expansion coefficients of Co(a) in the toy model have the same behavior. This suggests that non- 
perturbative corrections to T cuspid) and rcu°sp'(5') are parameterized by the same scale A^ defined 
in (I2.53p . Indeed we will show this in the next section by explicit calculation. 

We demonstrated in this section that nonperturbative corrections in the toy model follow 
unambiguously from the exact solution (I2.50p . In the next section, we will extend analysis to the 
cusp anomalous dimension and work out the strong coupling expansion of Tcuspls') / (2(7) analogous 
to (12391) . 
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3 Solving the quantization conditions 



Let us now solve the quantization conditions fl2.45p for the cusp anomalous dimension. The 
relation (12.451) involves two sets of functions. The functions Vo(t) and Vi(t) are given by the 
Whittaker function of 1st kind (see Eq. (ID. 21) ). At the same time, the functions fo(t) and fi(t) 
are defined in (I2.40p and they depend on the (infinite) set of expansion coefficients c{g) and 
c±{n,g). Having determined these coefficients from the quantization conditions (12.451) . we can 
then compute the cusp anomalous dimension for arbitrary coupling with a help of (I2.47p . 

We expect that at strong coupling the resulting expression for rcusplf?) will have the form (II. 2p . 
Examining (I2.47P we observe that the dependence on the coupling resides both in the expansion 
coefficients and in the functions U^{AiTg). The latter are given by the Whittaker functions of 2nd 
kind (see Eq. ( 1D.7P ) and, as such, they are given by Borel summable sign- alternating asymptotic 
series in 1/g. Therefore, nonperturbative corrections to the cusp anomalous dimension (12.470 
could only come from the coefficients c±{n,g) and c{g). 



3.1 Quantization conditions 

Let us replace fo{t) and /i(t) in (]2.45p by their explicit expressions ( I2.40p and rewrite the quan- 
tization conditions (I2.45P as 

VoiAngxi) + c{g)Vi{4:ngxe) = ^[c+{n, g)A+{n, xi) + c_{n, g)A_{n, xe)] , (3.1) 

re>l 

where Xg = i — ^ (with i = 0, ±1, ±2, . . .) and the notation was introduced for 

nVi{ATTgxe)Uo{A7rng) + XiVo{ATTgxe)U^{A7rng) 
A±{n,Xi) = . (3.2) 

The relation (13. ip provides an infinite system of linear equations for c±{g,n) and c{g). The 
coefficients in this system depend on Vo^i{ATTgXi) and U^i{ATrng) which are known functions 
defined in Appendix [Dl We would like to stress that the relation (13. ip holds for arbitrary g > 0. 

Let us show that the quantization conditions (13.11) lead to c{g) = for arbitrary coupling. To 
this end, we examine (13. ip for ^ 1. In this limit, for g = fixed we are allowed to replace the 
functions Vo{4:7rgxe) and VilAngXi) in both sides of (13. ip by their asymptotic behavior at infinity. 
Making use of flPTOj) and flDJ2D . we find for jx^l > 1 



, , VAAngxA ( -IQTrgXi + . . . , {xi < 0) 

'■<^'> " MSJJ? = I 1 + ..., (.,>o) 

where ellipses denote terms suppressed by powers of l/{gxi) and e"®'^^'!^^'. We divide both sides 
of (13. ip by Vi{4:7Tgx£) and observe that for X£ — * — oo the first term in the left-hand side of (13.10 
is subleading and can be safely neglected. In the similar manner, one has A±{n, xe) /VilAn gxe) = 
0(1/ xe) for fixed n in the right-hand side of (13.11) . Therefore, going to the limit — > — oo in 
both sides of (13.11) we get 

cig) = (3.4) 

for arbitrary g. We verify in Appendix El by explicit calculation that this relation indeed holds 
at weak coupling. 
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Arriving at ( I3.4p . we tacitly assumed that the sum over n in (I3.ip remains finite in the hmit 
xi — oo. Taking into account large n behavior of the functions U^{ATTng) and U^{ATrng) (see 
Eq. (1D.12P ). we obtain that this condition translates into the following condition for asymptotic 
behavior of the coefficients at large n 

c+(n, g) = o(n^/^) , c_(n, g) = o(n"^/^) . (3.5) 

These relations also ensure that the sum in the expression fl2.47p for the cusp anomalous dimension 
is convergent. 

3.2 Numerical solution 

To begin with, let us solve the infinite system of linear equations (13.11) numerically. In order to 
verify (13. 4p . we decided to do it in two steps: we first solve (13.11) for c±{n,g) assuming c{g) = 
and, then, repeat the same analysis by relaxing the condition (13. 4p and treating c{g) as unknown. 

For c{g) = 0, we truncate the infinite sums on the right-hand side of (13.11) at some large 
nmax and, then, use (13.11) for £ = 1 — n^^x, ■ ■ ■ ,'"'max to find numerical values of c±{n,g) with 
1 <n < n^i^^ for given coupling g. Substituting the resulting expressions for c±{n,g) into (12.471) 
we compute the cusp anomalous dimension. Taking the limit rimax — > cxd we expect to recover the 
exact result. Results of our analysis are summarized in two tables. Table 1 shows the dependence 
of the cusp anomalous dimension on the coupling constant. Table 2 shows the dependence of the 
cusp anomalous dimension on the truncation parameter nmax for fixed coupling. 

For c{g) arbitrary, we use (13.11) for i = —nmax, • • • , ""-max to find numerical values of c{g) and 
c±{n,g) with 1 < n < n^ax for given coupling g. In this manner, we compute Tcnsp{g)/{'2g) and 
c{g) and, then, compare them with the exact expressions corresponding to nmax oo. For the 
cusp anomalous dimension, our results for rcusp(fi')/(2(7) are in remarkable agreement with the 
exact expression. Namely, for nmax = 40 their difference equals 5.480 x 10~^ for g = 1 and it 
decreases down to 8.028 x 10"'' for g = 1.8. The reason why agreement is better compared to the 
c{g) = case (see Table 1) is that c{g) takes effectively into account a reminder of the sum in the 
right-hand side of (13. ip corresponding to n > nmax- The dependence of the obtained expression 
for c{g) on the truncation parameter nmax is shown in Table 3. We observe that, in agreement 
with (13.41) . c{g) vanishes as nmax — * oc. 



9 


0.1 


0.2 


0.4 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


1.8 


numer 


0.1976 


0.3616 


0.5843 


0.7096 


0.7825 


0.8276 


0.8576 


0.8787 


0.8944 


0.9065 


exact 


0.1939 


0.3584 


0.5821 


0.7080 


0.7813 


0.8267 


0.8568 


0.8781 


0.8938 


0.9059 



Table 1: Comparison of the numerical value of Tcnspid)/ i'^g) found from (13.11) and (I2.47p for 
'T-max = 40 with the exact one [25], HS] for different values of the coupling constant g. 

Our numerical analysis shows that the cusp anomalous dimension (12.470 can be determined 
from the quantization conditions (13.10 and (13.40 for arbitrary coupling g. In distinction with the 
toy model ( ]2.50p . the resulting expression for rcusp(5') does not admit a closed form representation. 
Still, as we will show in the next subsection, the quantization conditions (13. ip can be solved 
analytically for g ^ 1 leading to asymptotic expansion for the cusp anomalous dimension at 
strong coupling. 
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^max 


10 


20 


30 


40 


50 


60 


70 


oo 


numer 


0.8305 


0.8286 


0.8279 


0.8276 


0.8274 


0.8273 


0.8272 


0.8267 



Table 2: Dependence of Tcnspig) /{'^g) on the truncation parameter n„ 
The last column describes the exact result. 



for (7 = 1 and c{g) = 0. 



'^max 


10 


20 


30 


40 


50 


60 


70 


oo 


-c{9) 


0.0421 


0.0357 


0.0323 


0.0301 


0.0285 


0.0272 


0.0262 






Table 3: Dependence of c{g) on the truncation parameter 
quantization condition (13.11) . 



ioT g = 1 derived from the 



3.3 Strong coupling solution 

Let us divide both sides of (13.1 1) by Vo{4:TTgxe) and use (13.41) to get (for Xi 

. _Sr^ \nU(l{4:nng)r{xe) + U^{4:7mg)xe 

^ — c+(n, g) 

+ J^c_(n,^) 



i and ^ G Z) 



n — xi 

uUq {4:7mg)r{x£) + U{ {4:nng)x£ 



(3.6) 



n + xe 



where the function r{xg) was defined in (13. 3p . 

Let us now examine the large g asymptotics of the coefficient functions accompanying c±{n, g) 
in the right-hand side of (13.61) . The functions U^{4:7rng) and U^{A'Kng) admit asymptotic ex- 
pansion in 1/(7 given by flD.12p . For the function r{xi) the situation is different. As follows 
from its definition, Eqs. (13.31) and fID.lOp . large g expansion of r{x() runs in two parameters: 
perturbative 1/g and nonperturbative exponentially small parameter ~ (j,i/2g-27rg -^j-^iQj-^ -y^^g 
already encountered in the toy model, Eq. (12.531) . Moreover, we deduce from (13. 3p and fID.lOp 
that the leading nonperturbative correction to r{xi) scales as 



5r{x,)=0{A\^'-^\), 



e Z) 



(3.7) 



so that the power of A grows with i. We observe that O(A^) corrections are only present in 
r{xe) for £ = 0. Therefore, as far as the leading O(A^) correction to the solutions to (13.61) are 
concerned, we are allowed to neglect nonperturbative (A^— dependent) corrections to r(a;^) in the 
right-hand side of ( 13. 6p for i ^ and retain them for i = only. 

Since the coefficient functions in the linear equations (13.61) admit a double series expansion 
in powers of 1/g and A^, we expect that the same should be true for their solutions c±(n, g). Let 
us determine the first few terms of this expansion using the following ansatz: 



c±{n,g) = (Sngn) 



±1/4 



, , b±(n) 



a±[n) H : h . . 



+ 0(A^) , (3.8) 



where A^ is a nonperturbative parameter defined in (12.531) 

A2 = a(27r(7)i/2e-2-9^^i^ 



r(f)' 



(3.9) 
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and ellipses denote terms suppressed by powers of 1/g. Here the functions a±{n),h±{n), . . . are 
assumed to be (7— independent. We recall that the functions c±{n,g) have to verify the relation 
( 13.51) . This implies that the functions cL±{n),b±{n), . . . should vanish as n ^ 00. To determine 
them we substitute (13.81) into (13.61) and compare the coefficients in front of powers of 1/g and 
in both sides of (13.61) . 



3.3.1 Perturbative corrections 

Let us start with 'perturbative', A^— independent part of (13.81) and compute the functions a±(n) 
and b±{n). 

To determine a±{n), we substitute (13. 8p into (13. 6p . replace the functions UQ-^^{ATcgn) and r{xe) 
by their large g asymptotic expansion, Eqs. (ID. 121) and (13. 3p . respectively, neglect corrections in 
A^ and compare the leading 0{g^) terms in both sides of (13. 6p . In this way, we obtain from (13. 6p 
the following relations for a±(n) (with xi = i — ^) 

2x,T{l)J2^ = l, ii>l) (3.10) 

-2x,r(|)5^^ = l, (£<0) 
One can verify that the solutions to this system satisfying a±(ri) for n 00 have the form 



a_(ri) 



2r(n + l) 
n + l)l 
T(n + 



W = T./„ , iA-n2^iN ' (3-11) 



r(n + l)r2(i) 



4' 



2r(n + i)r2(: 



In the similar manner, we compare the subleading 0(1/ g) terms in both sides of (13.61) and find 
that the functions b±{n) satisfy the following relations (with x^ = i — j) 



2.,r(|)^i!±i!^ = --f -i^-i^i,>2. (^>i) (3.12) 

^-^ n — xp oZxp 64 32 

n>l 

-2x,r(f)VM^ = -^-^ + ^ln2, (£<0) 
^ ^^'^n + Xi 32xe 64 32 ' ^ - ; 

n>l 

where in the right-hand side we made use of (13. lip . Solutions to these relations are 

, . s . V31n2 3 \ , , 

b+{n) = -a+{n) ( + — ) , (3.13) 

b^{n) = a^{n) 

It is straightforward to extend analysis to subleading perturbative corrections to c±{n,g). 



31n2 




3 




+ 




4 


32r2 


31n2 




5 




+ 




4 


32n 
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Let us substitute (13. Sp into expression (\2A7\\ for the cusp anomalous dimension. Taking into 
account the identities ( ID. 121) we find the 'perturbative' contribution to rcusp(5') as 



r,,sp(^7) = 2g- ^(2 



Trn] 



n>l 



r( 



+r( 



a_M + + ... 



USngn 
3 



+ 



+ . . . 



(3.14) 



+ 0(A2). 



Attq J \ 1287cgn 

Replacing a±{n) and b±{n) by their expressions (13. lip and (13.130 . we find after some algebra 

31n2 K 



+ 0{l/g' 



(3.15) 



where K is the Catalan number. This relation is in agreement with the known result obtained 
both in TV = 4 SYM theory [251 126] and in string theory [29] . 

3.3.2 Nonperturbative corrections 

Let us now compute the leading O(A^) nonperturbative correction to the coefficients c±{n,g). 
According to (13. 8p . it is described by the functions a±{n) and (3±{n). To determine them from 
(13. 6p . we have to retain in r(x^) corrections proportional to A^. As was already explained, they 
only appear for i = 0. Combining together the relations (13. 3p . (ID.lOp and (1D.12P we find after 
some algebra 



4ng-- + 0{g-'] 



+ 0(A^ 



(3.16) 



Let us substitute this relation into (13. 6p and equate to zero the coefficient in front of A^ in the 
right-hand side of (13. 6p . This coefficient is given by series in 1/g and, examining the first two 
terms, we obtain the relations for the functions a±{n) and l3±{n). 

In this way, we find that the leading functions a±{n) satisfy the relations (with xe = i — j) 



2x,r(| 



-2x,r(|)5^ 



n — xp 

n>l 



0, 



TT 



n>l 



Je,o ■ 



(£>1) 
(£<0) 



where in the right-hand side we applied (13.110 . Solution to (13.170 satisfying a±{n) 
reads 



(3.17) 

as n — > oo 



a^{n) = , 

a^{n) = a_(n — 1) . 

with a_(n) defined in (13.110 . For subleading functions I3±{n) we have similar relations 

(3+{n) 1 



(3.18) 



2xpV{^ 



E 



n>l 



-2x,r(|)5^ 



n>l 



n — Xi 
n + Xi 



2 
1 



3tt 



8 16^2 



;i- 2 In 2)5^,0: 



ii>l) 
{i<0) 



(3.19) 
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In a close analogy with (]3.13p . the solutions to these relations can be written in terms of leading- 
order functions a±{n) defined in (13. lip 



= a_(n- 1) 



31n2 1 
—— + 



32n 



(3.20) 



It is straightforward to extend analysis and compute subleading O(A^) corrections to (13. 8p . 

The relation (13.80 supplemented with (13. lip . (13.130 . (13.180 and (I3.20p defines the solution to 
the quantization condition (13. 6p to leading order in both perturbative, 1/g, and nonperturbative, 
A^, expansion parameters. We are now ready to compute nonperturbative correction to the cusp 
anomalous dimension (I2.47p . Substituting (13. 8p into (12.470 we obtain 



5re,sp(^7) = -A'$^(27™)-4r(f)(', 

n>l L ^ 



0:4. H : h . . . 



+r(|)U_(n 



Airg 
Ang 



1 - 



+ . . . 



1 + 



1287r^n 
3 

1287r^n 



+ 



+ . . . 



(3.21) 



+ 0(A^). 



We replace a±{n) and P±{n) by their explicit expressions (13.180 and (I3.20p . evaluate the sums 
and find 



A2 

S'^cuspig) = 



TT 



3 - 6 In 2 
IGng 



+ 0{l/g' 



O(A^) 



(3.22) 



with A^ defined in (M . 

The relations (13.150 and (13.220 describe, correspondingly, perturbative and nonperturbative 
corrections to the cusp anomalous dimension. Let us define a new nonperturbative parameter 



cusp 



whose meaning will be clear in a moment 



m 



cusp 



4^2 



TTCr 



1 + ^^+0(1/9=) 

IbT^g 



+ 0(A^ 



(3.23) 



Then, the obtained expressions (13.150 and (I3.22p for the cusp anomalous dimension takes the 
form 



'2g 



31n2 

27T 



K 



8n^g 



+ 0{l/g' 



^"^cusp + 0{m] 



cusp/* 



(3.24) 



We recall that another nonperturbative parameter was already introduced in Sect. 2.1 as defining 
the mass gap mo(6) in the 0(6) model. We will show in the next section, that the two scales, 



m, 



cusp 



and mo (6), coincide to any order in 1/g. 



4 Mass scale 

The cusp anomalous dimension controls the leading logarithmic scaling behavior of the anomalous 
dimensions (11.10 in the double scaling limit L , N ^ 00 and j = L/ InN = fixed. The subleading 
corrections to this behavior are described by the scaling function e(j, g). At strong coupling, this 
function coincides with the energy density of the ground state of the bosonic 0(6) model (II. 3p . 
The mass gap in this model mo(6) is given by expression (12.90 which involves the functions T±{t) 
constructed in Section 2. 
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4.1 General expression 



Let us apply fl2.9p and compute the mass gap mo{6) at strong coupling. At large g the integral 
in (12.91) receives a dominant contribution from t ^ g. In order to evaluate (12.91) it is convenient 
to change the integration variable as t ^ Angit 



TT 



t + 



(4.1) 



where integration goes along the imaginary axis. We find from (I2.39P that T {An git) takes the 
form 



T{A7igit) = foiAngt)VoiA7igt) + MA7rgt)ViiAngt) , 



(4.2) 



where Vo^i{A7igt) are given by the Whittaker functions of first kind, Eq. (]D.2p . and fo^i{Angt) 
admit the following representation (see Eqs. (12.401) and (13.41) ) 



fo{Angt) = J2t 

n>l 

f,{Angt) = Y,r^ 



n>l 



n — t n + t 

n — t n + t 



- 1 



(4.3) 



Here the functions U^i {Airng) are expressed in terms of Whittaker functions of first kind, Eq. (]D.7p 
and the expansion coefficients c±{n,g) are solutions to the quantization conditions (I2.45p . 

Replacing T^Angit) in (14.10 by its expression (14.20 . we evaluate the t— integral and find after 
some algebra (see Appendix [E] for details) [38] 



"^0{6) 



16^2 



TT 



[foi-irgW^ing) + f,{-ng)U^ (ng)] 



(4.4) 



This relation can be further simplified with a help of the quantization conditions (I2.45p . For 
£ = 0, we obtain from (12.450 that /o(— 7r5')Vo(— vr^f) + fi{—TTg)Vi{—TTg) = 0. Together with 
the Wronskian relation for the Whittaker functions (ID. 80 this leads to the following remarkable 
relation for the mass gap 



"^0(6) 



16V2f\{-'Kg) 
7r2 Voi-irg) 



(4.5) 



It is instructive to compare this relation with similar relation (I2.46P for the cusp anomalous 
dimension. We observe that both quantities involve the same function fi{Angt) but evaluated 
for different values of its argument, that is t = —1/4 for the mass gap and t = for the 
cusp anomalous dimension. As a consequence, there are no reasons to expect that the two 
functions, m{g) and Tcnsp{g), could be related to each other in a simple way. Nevertheless, we 
will demonstrate in this subsection, that ^o(6) determines the leading nonperturbative correction 
to Fcuspls') at strong coupling. 
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4.2 Strong coupling expansion 

Let us now determine the strong coupling expansion of the functions (14.31) . We replace coefficients 
c±{n,g) in (14.31) by their expression (13. 8p and take into account the obtained results for the 
functions a±,b±,..., Eqs. (I3.1ip . (13.131) . (13.181) and (I3.20p . In addition, we replace in (14.31) the 
functions U^i{ATTng) by their strong coupling expansion ( ID. 121) . We recall that the coefficients 



c±in,g) admit the double series expansion ( 13.81) in powers of 1/g and A ~ e" 
a consequence, the functions folAngt) and fi{4TTgt) have the form 

UiArrgt) = fr^'^^Qt) + ^U^T^Qt) , (n = 0, 1) 



-2ng 



Eq. (EH). As 



(4.6) 



where /^^^' is given by asymptotic (non-Borel summable) series in 1/g and 5/„ takes into account 
nonperturbative corrections in A^. 

Evaluating sums in the right-hand side of (14. 3p we find that f 0(4:71 gt) and fi{4:Tigt) can be 
expressed in terms of two sums involving functions a±{n) defined in ( 13. lip 



2r( 



2r( 



a_|_(n) 



^' ^ t-n 

n>l 



E 

n>l 



a_{n) 
t + n 



r(|)r(i^ 
t) 



r(! 



T(i)r(i + t) 
r(i + t) 



(4.7) 



Going through calculation of (14. 3p . we find after some algebra that perturbative corrections to 
fo{4:7rgt) and fi{4Trgt) are given by linear combinations of the ratios of Euler gamma-functions 



r(r 



r(i-t) 



+ 



1 



Aug 
1 

Ang 
1 



3 In 2 



+ 



i\ r(|)r(i-t) r(|)r(i + t) 
it) r(|-t) 8tT{l + t) 



+ 0(^7" 



r(i)r(i + t) r(f)r(i-t) 
4tr(i + t) 
r(i)r(i + t) 



4tr(f-t) 



(4.8) 



(ATcg)^ [ Atra + t) 



1 

4t 



31n2 



r(|)r(i 



t) 



AtTU-t) 



1 31n2 

At ^ A 



+ 0{g- 



Notice that fi{t) is suppressed by factor l/{A7ig) compared to fo{t)- In the similar manner, we 
compute nonperturbative corrections to (14. 6 p 



6foiAngt) = A^ 
5f,{A7Tgt) = A' 



A-ng 



r(|)r(i 



t) 



r(r 



T{l + t) 



2r(f-t) 



2r(f + t) 



+ 0{g 



-2\ 



f 1 r(|)r(i + t) 
\A7Tg r(f + t) 



(4.9) 



+ 



{Airgy 



T(|)r(i 



t) 



^) f ^ 3, ^ 1 
^— h - In 2 

rf + t \%t A A 



r(f)r(i 



+ 0{g 



-3n 



where ellipses denote O(A^) terms. 
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Substituting (I4.8P and (I4.9P into (I4.2p we obtain the strong coupling expansion of the function 
T^inigt). To verify the obtained expressions, we apply fl2.46p to calculate the cusp anomalous 
dimension 

rcusp(^7)=2^7-4(7/r'(0)-4(?5/i(0). 
Replacing f[^^\0) and 5/i(0) by their expressions, Eqs. (14. Sp and fl4.9p . we obtain 



(4.10) 



cusp 



3 In 2 



K 



+ ... 



A2 

TT 



1 + 



3 - 61n2 
16ng 



+ ... 



(4.11) 



in a perfect agreement with ( 13.15^ and ( ]3.22p . respectively. 

Let us obtain the strong coupling expansion of the mass gap (14. 5p . We replace Vo(— vr^f) by 
its asymptotic series, Eqs. ( 1D.14P and ( ]D.12p . and take into account (14. 8 p and (14. 9p to get 



"^0(6) 



V2_ 

r(4) 



(27r(7) 



1/4, 



-Tvg 



3-61n2 -63 + 108 ln2 - 108(ln2)2 + 16K 

1 + + + ... 
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A2 



Sng 



2048(7r^)2 
15-61n2 
327r7 



+ . . . 



+ 0(A^)L (4.12) 



Here, in order to determine 0{1/ g"^) and 0{A'^/g'^) terms inside the curly brackets, we computed 
in addition the sub leading 0{g~^) corrections to /{^^' and 6fi in Eqs. (14. 8 p and (14. 9p . respectively. 
The leading 0(1/5^) correction to mo(6) (the second term inside the first square bracket in the 
r.h.s. of (I4.12P ) is in agreement with both analytical [35l [38] and numerical calculations |36j . 

We are now ready to clarify the origin of the 'substitution rule' (I2.66P that establishes the 
relation between the cusp anomalous dimension in the toy model and the exact solution. To 
this end, we compare the expressions for the functions fniAngt) given by (14. 6p . (14. 8 p and (14.90 
with those in the toy model, Eqs. (I2.42p and (l2.48p FI It is straightforward to verify that upon 
the substitution (12.660 and (I2.64p the two set of functions coincide up to an overall t— dependent 
factoi0 



fni^ngt) 



r(|)r(i 



(n 



0,1) 



(4.13) 



Since the cusp anomalous dimension (I2.46P is determined by the /i— function evaluated at t = 0, 
the additional factor does not affect its value. 



4.3 Nonperturbative corrections to the cusp anomalous dimension 

The relation (14.120 defines strong coupling corrections to the mass gap. In a close analogy 
with the cusp anomalous dimension (14. lip , it runs in two parameters: perturbative 1/g and 
nonperturbative A^. We would like to stress that the separation of the corrections to mo(6) into 
perturbative and nonperturbative ones is ambiguous since the 'perturbative' series inside the 
square brackets in the right-hand side of ( 14.120 is non-Borel summable and, therefore, it suffers 
from Borel ambiguity. It is only the sum of perturbative and nonperturbative corrections that 

^It worth mentioning that the functions /o*°'*'^ and /l*"^-* in the toy model are, in fact, i— independent. 

^Roughly speaking, this substitution simpUfies the comphcated structure of poles and zeros of the exact 
solution, Eqs. (|4.8p and (14. 9p . encoded in the ratio of the gamma- functions to match simple analytical properties 
of the same functions in the toy model (compare (|2.13p and (|2.16p ). 



28 



is a unambiguously defined function of tlie coupling constant. In distinction with the mass scale 
mo(6) ; the definition (12.531) of the nonperturbative scale involves a complex parameter a whose 
value depends on the prescription employed to regularize singularities of the 'perturbative' series. 

To illustrate the underlying mechanism of cancellation of Borel ambiguity inside mo(6), let 
us examine the expression for the mass gap (14. 5 p in the toy model. As was already explained in 
Sect. 2.8, the toy model captures the main features of the exact solution at strong coupling and, 
at the same time, it allows us to obtain expressions for various quantities in a closed analytical 
form. The mass gap in the toy model is given by the relation (14. 5 p with fi{—7fg) replaced with 
fi''^\-TTg) defined in (E^Sl) and fICTD . In this way, we obtain 

16v/2/^)(-vr,) 16V2 1 

Here Vi(— vr^f) is an entire function of the coupling constant (see Eq. (12.411) ). Its large g asymptotic 
expansion can be easily deduced from flD.16p and it involves the nonperturbative parameter A^. 
Making use of (I2.52p we obtain from (14.140 



^toy = ^(2vr^7)^/^e--^| 



1 23 

+ 



327r^ 2048(7r^)2 

A2 r 11 

1 h 

Sug 32ng 



+ 0(A^) , (4.15) 



where ellipses denote terms with higher power of 1/g. By construction, mtoy is a unambiguous 
function of the coupling constant whereas the asymptotic series inside the square brackets are 
non-Borel summable. It is easy to verify that 'perturbative' corrections to w^toy described by 
the asymptotic series C2{a) given by (I2.6ip . Together with (I2.59P this allows us to identify the 
leading nonperturbative correction to (I2.59P in the toy model as 

A 2 2 

= --C2(«) + O(A^) = -^aml^ + 0{mtj , (4.16) 

with A^ given by fl3^ . 

Comparing the relations (I4.12p and (I4.15P we observe that mo(6) and mtoy have the same 
leading asymptotics while subleading 1/g corrections to the two scales have different transcen- 
dentality. Namely, the perturbative coefficients in mtoy are rational numbers while for mo(6) their 
transcendentality increases with order in 1/g. We recall that we already encountered the same 
property for the cusp anomalous dimension, Eqs. (I2.59P and (I2.63p . There, we have observed 
that the two expressions (I2.59P and (12.630 coincide upon the substitution (12.640 . Performing the 
same substitution in (I4.12p we find that, remarkably enough, the two expressions for the mass 
gap indeed coincide up to an overall normalization factor 

Eq.lMl TT 

The expressions for the cusp anomalous dimension (14. lip and for the mass scale (14. 12p can 
be further simplified if one redefines the coupling constant as 

31n2 

g' = g-ci, ci = -— , 4.18 

4:71 
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and re-expands both quantities in 1/g'. As was observed in |25|, such redefinition allows one to 
eliminate 'In 2' terms in perturbative expansion of the cusp anomalous dimension. Repeating the 
same analysis for fl4.1ip we find that the same is also true for nonperturbative corrections 



Tcusp {g + ci) =2g 



1 - 



K 



+ ... 



A2 



2v^: 



TT 



1 + 



IQixg 



(4.19) 



with A^ defined in (13. 9p . In the similar manner, the expression for the mass scale (14.121) takes 
the form 



[mo(&){g + ci 



2A2 



TTCr 



3 16K - 54 



IQ'Kg 512(7r^)2 



+ 



+ 0(A^) 



(4.20) 



Comparing the relations (14.19P and (14.201) we immediately recognize that, within an accuracy of 
the obtained expressions, nonperturbative ©(A^) correction to the cusp anomalous dimension is 
given by 



o(6) 



5r, 



cr 



cusp 



4^2 



o(6) 



+ 0(m^(6)). 



(4.21) 



It worth mentioning that, upon identification of the scales (I4.17p . this relation coincides with 

We will show in the next subsection that the relation (I4.2ip holds at strong coupling to all 
orders in 1/g. 



4.4 Relation between cusp anomalous dimension and mass gap 

We demonstrated that the strong coupling expansion of the cusp anomalous dimension has the 
form (I3.24P with the leading nonperturbative correction given to the first few orders in 1/g 
expansion by the mass scale of the 0(6) model, m1^^^ = m^f^^y Let us show that this relation is 
in fact exact at strong coupling. 

According to (I4.10p . the leading nonperturbative correction to the cusp anomalous dimension 
is given by 

5rcusp = -4(7 5/i(0), (4.22) 

with 5/1(0) denoting 0{A'^) correction to the function fi(t = 0), Eq. (14. 6p . We recall that this 
function verifies the quantization conditions (12.450 . As was explained in Section 3.3, the leading 
0(A2) corrections to solutions of (I2.45P originate from subleading, exponentially suppressed terms 
in the strong coupling expansion of the functions Vo(— vrgf) and Vi{—'Kg) that we shall denote as 
SVo^—ng) and 6Vi{—TTg), respectively. Using the identities flD.14p and flD.lOp . we find 

SVo{~7rg) = a'^e--3u-^^g^^ (4.23) 
SVi{-7rg) = a^e--^Uring), 
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where the functions UQ^ng) and [/{{ng) are defined in (1D.7I) . Then, we spht the functions fo(t) 
and fi{t) entering the quantization conditions (12.451) into perturbative and nonperturbative parts 
according to (14. 61) and compare exponentially small terms in both sides of fl2.45l) to get 

SfoikWoik) + 6f\{te)Vi{t,) = -m'6e,o , (4.24) 

where ti = Angi^i — jj and the notation was introduced for 

m' = fo{-7ig)SVo{-ng) + f,(-ng)SVi{-ng) . (4.25) 

Taking into account the relations (I4.23P and comparing the resulting expression for m' with (14.41) 
we find that 

,0" 

m = -— mo{6) , (4.26) 

with mo(6) being the mass scale (14. 4p . 

To compute nonperturbative O(A^) correction to the cusp anomalous dimension, we have to 
solve the system of relations (I4.24p . determine the function Sfi(t) and, then, apply (I4.22p . We 
will show in this subsection that the result reads 

(5/1(0) = — ^m'mo(6) = o— m^(6) , (4.27) 

to all orders in strong coupling expansion. Together with (I4.22p this leads to the desired expres- 
sion (I4.2ip for leading nonperturbative correction to the cusp anomalous dimension. 
To begin with, let us introduce a new function analogous to (I2.39p 

bV{it) = 6fo{t)Vo{t) + 6Mt)Vi{t) . (4.28) 

Here Sfo{t) and Sfi{t) are given by the same expressions as before, Eq. (I2.40p . with the only 
difference that the coefficients c±{n,g) are replaced in (I2.40p by their leading nonperturbative 
correction Sc±{n,g) = O(A^) and the relation (13. 4p is taken into account. This implies that 
various relations for T{it) can be immediately translated into those for the function 6r{it). In 
particular, for t = we find from (12.400 that 5/o(0) = for arbitrary coupling, leading to 

5r(0) = 25/i(0) (4.29) 

In addition, we recall that, for arbitrary c±{n,g), the function (I2.39P satisfies the inhomogeneous 
integral equation (12. 7p . In other words, the c-i-(?t,, (7)— dependent terms in the expression for the 
function T{it) are zero modes for the integral equation (12. 7p . Since the function (14.280 is just 
given by the sum of such terms, it automatically satisfies the homogenous equation 



dt 



e**" (5r_(t) -e~**"5r+(t) 



0, (-1<M<1), (4.30) 



where 6T(J:) = 6T+{t) + i6T_{t) and 6T±{-t) = ±6T±{t). 

As before, in order to construct the solution to (I4.30p . we have to specify additional conditions 
for Sr{t). Since the substitution c±{n,g) — * Sc±{n,g) does not affect analytical properties of the 
functions (12.400 . the function (I4.28P shares with T{it) an infinite set of simple poles located at 
the same position (I2.14p 

5^) ~ ^-1^ , {i G Z/0) . (4.31) 
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In addition, we deduce from ( I4.24p that it also satisfies the relation (with 

6T{4:TTigx£) = —m'6efi , (£ G Z) , 



(4.32) 



and, therefore, has an infinite number of zeros. An important difference with T{it) is that ST (it) 
does not vanish at t = —irg and its value is fixed by the parameter m' defined in fl4.26p . 

Having in mind similarity between the functions T{it) and 6r{it) we follow fl2.13p and define 
a new function 

= ^ ""(^/'^^ 5^) . (4.33) 

^ ' V2sin(t/4(7 + 7r/4) ^ ' ^ ' 

As before, the poles and zeros of T{it) are compensated by the ratio of sinus functions. However, 
in distinction with '~f{it) and in virtue of 5r{—7rig) = —m', the function S'y{it) has a single pole 
aX t = —TTQ with the residue equal to 2gm'. For t we find from fl4.33p that 6'~f{it) vanishes as 



t 



t 



(4.34) 



6^{zt) = -5r(0) + 0{t') = -Sf\{0) + 0{t') , 

where in the second relation we applied (14.291) . It is convenient to split the function 6'y{t) into the 
sum of two terms of a definite parity, S'jit) = S'y+(t) + i6j_(t) with S'y±{—t) = ±5'y±(t). Then, 
combining together (I4.30p and (I4.33P we obtain that the functions S'y±{t) satisfy the infinite 
system of homogenous equations (for n ^ 1) 



dt 



dt 



J2n{t) 



5^.{t) , (57+ (t) 



1 _ Q-t/2g 



+ 



Qt/2g _l 



-t/2g Qt/2g _l 



= 0, 

0. (4.35) 
and have a simple pole 



By construction, the solution to this system 6'y{t) should vanish at t 
at t = —ing. 

As was already explained, the functions 5T±{t) satisfy the same integral equation f l4.30p as 
the function T±{t) up to inhomogeneous term in the right-hand side of (12. 7p . Therefore, it should 
not be surprising that the system (I4.35P coincides with the relations (12.30 after one neglects the 
inhomogeneous term in the right-hand side of (12. 3p . As we show in Appendix [Cl this fact allows 
us to derive Wronskian like relations between the functions 57(t) and 7(t). These relations turn 
out to be powerful enough to determine the small t asymptotics of the function 57(t) at small t 
in terms of 7(t), or equivalently r(t). In this way we obtain (see Appendix O for more detail) 



5'~^{it) = —m't 



-wg 



ir'^g 



V2 



IT 



-irg 



Re 



dt' 



(4.36) 



'0 t' + iirg 

Comparing this relation with (12.90 we realize that the expression inside the square brackets is 
proportional to the mass scale mo(6) leading to 

tV2 



6'^ {it) 



-mmo(6)- 



8^7 



0{i' 



(4.37) 



Matching this relation into (14.340 . we obtain the desired expression for 5/i(0), Eq. (I4.27p . Then, 
we substitute it into (14.220 and compute the leading nonperturbative correction to the cusp 
anomalous dimension, Eq. (14.210 . leading to 



"^cusp(fi') = mo(Q){g) ■ 



(4.38) 
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Thus, we demonstrated in this section that nonperturbative, exponentially small corrections to 
the cusp anomalous dimensions at strong coupling are determined to all orders in 1/g by the 
mass gap of the two-dimensional bosonic 0(6) model embedded into AdSs x sigma-model. 

5 Conclusions 

In this paper, we have studied anomalous dimensions of Wilson operators in the SL{2) sector 
of planar Af = 4 SYM theory in the double scaling limit when Lorentz spin of the operators 
grows exponentially with their twist. In this limit, the asymptotic behavior of the anomalous 
dimensions is determined by the cusp anomalous dimension Tcnspig) and the scaling function 
^{QiJ)- We found that at strong coupling both functions receive exponentially small corrections 
which are parameterized by the same nonperturbative scale. It is remarkable that this scale 
appears on both sides of the AdS/CFT correspondence. In string theory it emerges as the mass 
gap of the two-dimensional bosonic 0(6) sigma model which describes the effective dynamics of 
massless excitations for folded spinning string in the AdSs x sigma model [6J. 

The dependence on rcusp(5') a-^id ^{g,j) on the coupling constant is governed by integral 
BES/FRS equations which follow from the conjectured all-loop integrability of the dilatation 
operator of A/" = 4 model. At weak coupling, their solutions agree with results of explicit pertur- 
bative calculations. At strong coupling, systematic expansion of the cusp anomalous dimension 
in powers of 1/g was derived in p5|. In agreement with the AdS/CFT correspondence, the 
first few terms of this expansion coincide with the energy of the semiclassically quantized folded 
spinning strings. However, the expansion coefficients grow factorially at higher orders and, as a 
consequence, 'perturbative' 1/g expansion of the cusp anomalous dimension suffers from Borel 
singularities which induce exponentially small corrections to rcusp(fl')- To identify such nonper- 
turbative corrections, we revisited the BES equation and constructed the exact solution for the 
cusp anomalous dimension valid for arbitrary coupling constant. 

At strong coupling, we found that the obtained expression for rcusp(fi') depends on a new 
scale mcusp(fl') which is exponentially small as g ^ oo. Nonperturbative corrections to rcusp(5') at 
strong coupling run in even powers of this scale and the coefficients of this expansion depend on 
the prescription employed to regularize Borel singularities in perturbative 1/g series. It is only 
the sum of perturbative and nonperturbative contributions which is independent on the choice 
of the prescription. For the scaling function e{g,j), the defining integral FRS equation can be 
brought to the form of the thermodynamical Bethe ansatz equations for the energy density of 
the ground state of the 0(6) model. As a consequence, nonperturbative contribution to e{g,j) 
at strong coupling is described by the mass scale of this model Tno{e){g). We have shown that 
the two scales coincide, rricuspig) = ""^0(6) (5'); and, therefore, nonperturbative contributions to 
rcusp(5') and e{g,j) are governed by the same scale mo(6)(fl')- 

This result agrees with the proposal by Alday-Maldacena that, in string theory, the leading 
nonperturbative corrections to the cusp anomalous dimension coincide with those to the vacuum 
energy density of two-dimensional bosonic 0(6) model embedded into the AdSs x sigma- 
model. These models have different properties: the former model has asymptotic freedom at 
short distances and develops mass gap in the infrared while the latter model is conformal. The 
0(6) model only describes an effective dynamics of massless modes of AdSs x and the mass 
of massive excitations /i ~ 1 defines a ultraviolet (UV) cut-off for this model. The coupling 
constants in the two models are related to each other as ^^(/i) = l/{2g). The vacuum energy 
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density in the 0(6) model and, more generally in the 0{n) model is a ultraviolet divergent 
quantity. It also depends on the mass scale of model and has the following form 

evac = f^'^^if) + i^ml^n) + 0{m^o{n)/f^^) ■ (5-1) 
Here //^ is a UV cut-off, e(^^) stands for perturbative series in g'^ and the mass gap ^o(n) 

ml(n) = cf^'e'^ r^P-l^l [1 + 0{f)\ , (5.2) 

where /3o and /3i are the beta-function coefficients for the 0(?7.) model and the normalization factor 
c ensures independence of mo(n) on the renormalization scheme. For n = 6 the relation (15.21) 
coincides with (11.31) and the expression for the vacuum energy density (15.11) should be compared 
with 

The two terms in the right-hand side of (15.11) describe perturbative and nonperturbative 
corrections to evac- For n — > oo each of them is well-defined separately and can be computed 
exactly [IHl US]. For n finite, including = 6, the function e(^^) is given in a generic renor- 
malization scheme by a non-Borel summable series and, therefore, is not well-defined. In a close 
analogy with (II. 2p . the coefficient k in front of '^o{n) right-hand side of (15.11) depends 

on the regularization of Borel singularities in perturbative series for e((7^). Notice that evac is 
related to the vacuum expectation value of the trace of the tensor energy-momentum in the 
two-dimensional 0(n) sigma model [19]. The AdS/CFT correspondence implies that for n = 6 
the same quantity defines nonperturbative correction to the cusp anomalous dimension (II. 2p . It 
would be interesting to obtain its dual representation (if any) in terms of certain operators in 
four- dimensional A/" = 4 SYM theory. Finally, one may wonder whether it is possible to iden- 
tify a restricted class of Feynman diagrams in A/" = 4 theory whose resummation could produce 
contribution to the cusp anomalous dimension exponentially small as (7 — * cxo. As a relevant 
example, we would like to mention that exponentially suppressed corrections were obtained in 
Ref. [50] from exact resummation of ladder diagrams in four-dimensional massless g(^^ theory. 
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A Weak coupling expansion 

In this Appendix, we work out the first few terms of the weak coupling expansion of the coefficient 
c{^g) entering (I2.47P and show that they vanish in agreement with (13.40 . To this end, we will not 
attempt at solving the quantization conditions (I2.45P at weak coupling but will use instead the 
fact that the BES equation can be solved by iteration of the inhomogeneous term. 

The system of integral equation (12.31) can be easily solved at weak coupling by looking for 
its solutions 7±(t) in the form of the Bessel series (12.101) and expanding the coefficients and 
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72fc-i in powers of the coupling constant. For g —>■ it follows from (I2.3P and from orthogonality 
conditions for the Bessel functions, that 7-(t) = Ji{t) + . . . and 7+(t) = + . . . with ellipses 
denoting subleading terms. To determine such terms it is convenient to change the integration 
variable in (12.31) as t — *■ tg. Then, taking into account the relations Jk{—gt) = {—l)''Jk{gt) we 
observe that the resulting equations are invariant under substitution g —g provided that the 
functions 'y±{gt) change sign under this transformation. Since 7±(— t) = ±7±(t), this implies 
that the coefficients ■j2n-i{g) and 'j2n{g) entering fl2.10p have a definite parity as functions of the 
coupling constant 

l2n-l{-g) = l2n~l{g) , 72n(-^) = -72n(^) , (A.l) 

and, therefore, their weak coupling expansion runs in even and odd powers of g, respectively. 
Expanding both sides of (12. 3p at weak coupling and comparing the coefficients in front of powers 
of g we find 

71 = 9 + 9 - + 4Cs 9 + O(g^) , 

' 2 6 ^ 90 ^ V 630 

72 = C39' - (yC3 + + (^Ca + ^C5 + 105Cr)^7^ + 0(/) , 

73 = -^9' + ^^9' + 0{g') , 74 = - (yC5 + 21Cr)/ + 0{g') , 

75 = + , 76 = C79' + 0{g') . (A.2) 

945 

We verify with a help of ( 12.12^ that the expression for the cusp anomalous dimension 

= ^9'lM = 4/ - ^9' + - + 32(1 jg' + 0(^^°) (A.3) 

agrees with the known four loop result in planar A/" = 4 SYM theory [20] . 

In our approach, the cusp anomalous dimension is given for arbitrary value of the coupling 
constant by the expression (12.471) which involves the functions c{g) and c±{n,g). According 
to (12.291) . the latter functions are related to the functions 7(t) = 7+(t) + il~{t) evaluated at 
t = Airign 

c+{n,g) = -4^e-^"3" [j+{4nign) + ij^{4mgn)] , 

c-{n,g)= Age"^'^^"' [y+{4TTign) - i^^{ATTign)] . (A.4) 

At strong coupling, we determined c±{n, g) by solving the quantization conditions (13. 6p . At weak 
coupling, we can compute c±{n,g) from flA.4p by replacing 7±(t) with their Bessel series (I2.10p 
and making use of the obtained expressions for the expansion coefficients (IA.2p . 

The remaining function c{g) can be found from comparison of two different representations 
for the cusp anomalous dimension, Eqs. (12.470 and (I2.12p . 

c{g) = + 2(?7i(^7) + Uoi^nng) + c+{n,g) U+{47ing)] . (A.5) 

n>l 
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Taking into account the relations (1A.4I) and (]2.10p we find 

c{g) = -\ + 2^71 (^7) - - i)i2k-i{g) f2k~i{g) + {mi2k{g)f2k{g)\ , (A.6) 

k>l 

where the coefficients 7^ are given by flA.2p and the notation was introduced for the functions 
fk{g) = SgY,M(^^9n) - {-lfU^{ATign)\h{A7:gn)e-^^^^ . (A.7) 

n>l 

Here Ik{x) is the modified Bessel function |17] and the functions Uq{x) are defined in flD.7p . At 
weak couphng, the sum over n can be evaluated with a help of the Euler-Maclaurin summation 
formula. Going through lengthy calculation we find 

/i = 1 - 2^ + + 2C,g' - %-g' - 23C,g' + + ^Cr/ + 0(/) , 

6 1U8 4 

/2 = + 2Csg' - ^g' + 0{g') , h = ^ + 0{g') , 

/4 = -^ + 0(/), h = l + 0{g'), /6 = -^ + 0(/). (A.8) 

In this way, we obtain from flA.Gp 

c{g) = -\ + ifi + 2^7)71 + 2/272 - 3/373 - 4/474 + 5/575 + 6/676 + . . . = . (A.9) 

Thus, in agreement with ( 13. 4p . the function c{g) vanishes at weak coupling. As was shown in 
Sect. 3.1, the relation c{g) = holds for arbitrary couphng. 



B Constructing general solution 

By construction, the function T(t) = T^(t) + iT^(t) defined as the exact solution to the integral 
equation (12. 7p is given by the Fourier integral 

/oo ^ 

dke-'^'T{k), (B.l) 
-00 

with the function r(A;) having different form for k"^ < 1 and k"^ > 1: 

• For —00 < k < —1: 

f{k)=J2c-{n,g)e-'--'^-'-'K (B.2) 

n>l 

• For 1 < k < 00: 

Tik)=J2c+in,g)e-'--^^'-'\ (B.3) 

n>l 
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For -1< A: < 1: 



m 



' -K \l-k 



1 + 



c{9) 



1 

k^2 



— 1 /"OO 

+ 

OO J 1 



dpT{p) f p — 1 



p 



p+l 



1/4 



(B.4) 



where T{p) inside the integral is replaced by (\B.2\i and fIB.Sp . 



Let us split the integral in fIB.ip into three terms as in fl2.38p and evaluate them one after another. 
Integration over fc^ > 1 can be done immediately while the integral over — 1 < < 1 can be 
expressed in terms of special functions 



n>l 

+ ^c_{n,g) 



Aimg + it 



V+i-it^Anng) 



n>l 



+ V_{it,ATxng) 



Vo{-it) - c{g)V,{-it) 



Aimg — it 

where the notation was introduced for the functions (with ri = 0, 1) 

1 r dpe-y^P~^Ul + kp-l\^^'^ 



(B.5) 



Vn{x) 



. dke , 

V27T J -I J I p-k 

V2 dke""^ fl + k^^'^ 



l-kp+l 



vr 7„i (fc + 1)" \l-k 



dp e-2^(P-i) (p + l 



Tl/4 



(B.6) 



The reason why we also introduced U^{y) is that the functions V±{x, y) can be further simplified 
with a help of master identities (we shall return to them in a moment) 

{x + y)V4x, y) = xVo{x)U{{y) + yVi{x)U^{y) - e"" , 

(x - y)V^{x, y) = xV,{x)Ut{y) + yV^{x)U^{y) - e^ . (B.7) 

Combining together ( 1B.7I) and ( IB.SP we arrive at the following expression for the function T{it) 

V{it) = ~V,{t) ~ c{g)V^{t) 

' A-nngVi{t)U^ {Aimg) + tVo{t)Ut {"^img)' 



+ ^c+{n,g) 

n>l 

+ ^c_(n,^) 



n>l 



Aiing — t 

ATrngVi{t)UQ {Airng) + tVo{t)U{ (Anng) 
Anng + 1 



(B.8) 



which leads to ^M)- 

We show in Appendix [D] that the functions Vo,i(t) and U^i{ATcng) can be expressed in terms 
of Whittaker functions of the first and second kind, respectively. As follows from their integral 
representation, Vo{t) and Vi(t) are holomorphic functions of t. As a result, T{it) is a meromorphic 
function of t with (an infinite) set of poles located at t = ^Anng with n positive integer. 
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Let us now prove the master identities (1B.7I) . We start with the second relation in (IB .71) and 



make use of (IB. 61) to rewrite the expression in the left-hand side of (IB. 71) as 

POO 

(x - y)V+{x, y) = (x - y) ds Voix + s)U^{y + s) e"^"^ . (B.9) 

^0 

Let us introduce two auxiliary functions 

zi{x) = Vi{x), zi{x) + z[{x) = Vo{x) , 

Z2{x) = e-^Ut{x) , Z2{x) + z'^ix) = -e-"?7o+(x) , (B.IO) 

with Vn{x) and f/^(a;) given by (IB. 60 . They satisfy the second-order differential equation 

^(xzlix)) = {x-^)z,{x). (B.ll) 

Applying this relation it is straightforward to verify the following identity 

-{x - y) [zi{x + s) + z[{x + s)] [z2iy + s) + z'^iy + s) 



= j;{(y + ^)[^2(y + s) + z'^iy + s)]zi{x + s)} 

- ^{(x + s)[zi{x + s) + z[{x + s)]z2iy + s)}. (B.12) 

It is easy to see that the expression in the left-hand side coincides with the integrand in (IB. 91) . 
Therefore, integrating both sides of ( IB. 121) over < s < oo, we obtain 



=oo 
=0 



{x - y)V+{x, y) = - e-' [{x + s)Vo{x + s)Ut{y + s) + {y + s)Vi{x + s)U+{y + s)] 

= -e- + xVo{x)Ut{y) + yVi{x)U+{y) , (B.13) 

where in the second relation we took into account the asymptotic behavior of the functions (IB. 61) 
(see Eqs. (IDAOD and (ID12D ). ~ s'^/^ and U+{s) ~ s""^/^ as s ^ oo. 

The derivation of the first relation in (1B.7P goes along the same lines. 

C Wronskian like relations 

In this Appendix we present a detailed derivation of the relation (14.361) which determines the 
small t expansion of the function 57 (t). This function satisfies the infinite system of integral 
equations (I4.35p . In addition, it should vanish at the origin, t = and have a simple pole at 
t = —iirg with the residue 2igm' (see Eq. (I4.33P ). To fulfill these requirements, we split 5'^{it) 
into the sum of two functions 

2m' f 

57(^t)=7(^t)-— — — , (C.l) 
TT t + Tcg 

where, by the construction, ^{it) is an entire function vanishing at t = and its Fourier transform 
has a support on the interval [—1, 1]. Similarly to (12.20 . we decompose 6'y{t) and 7(t) into the 
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sum of two functions with a definite parity 

57+(t)=7+(t) 
57_(t) =7_(t) + 



2m' 



2gm't 



t2 + 7r2^2 



(C.2) 



Then, we substitute these relations into fl4.35p and obtain the system of inhomogeneous integral 
equations for the functions 7±(t) 



T 
T 



J2n{t) 



+ 



1 _ e-*/2g Qt/2g _i 



1 



-t/2g 



with inhomogeneous terms given by 

2m' rdtJ2n-i{t) 



h 



2n-l 



h 



Jo + {'^aY 

2m' dtJ2n{t) 



ot/2g 



t 



2n 



et/{2g) 
Qt/{2g) _l 



+ 



= h2n-l{9) 
-- h2n{g) , 

1 _ Q~t/{2g) 
t 

1 _ e-</(29) 



(C.3) 



(C.4) 



Comparing these relations with flC.3P we observe that they only differ by the form of inhomoge- 
neous terms and can be obtained one from another through the substitution 



In a close analogy with (12.101) . we look for solution to ( 1C.3|) in the form of Bessel series 

7-(t) = 2^ (2r2 - l)J2n-l{t)l2n-l{9) , 

7+(t) = 2 J] (2n) J2„(t)72n(^?) . 
For small t we have 7_(t) = f^i + 0(t2) and 7+(t) = 0(t2). Then it follows from flUTT]) 



(C.5) 



(C.6) 



2m! 

-n'^g 



(C.7) 



so that the leading asymptotics is controlled by the coefficient 71. 

Let us multiply both sides of the first relation in (1C.3P by (2n — l)72n-i and sum both sides 
over n > 1 with a help of fl2.10p . In the similar manner, we multiply the second relation in (1C.3P 
by (2n)72n and follow the same steps. Then, we subtract the second relation from the first one 
and obtain 



7-(t)7-W-7+W7+W , 7-W7+W+7+W7-W 



+ 



1 _ Q-tl2g Qt/2g _l 

= 2 ^ [{2n - l)72n-l/i2n-l " (2n)72„/l2n] • 



n>l 
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We notice that the expression in the left-hand side of this relation is invariant under exchange 
7±(^) ^ 7±(^)- Therefore, the right-hand side should be also invariant under (jC.511 leading to 



7i 



2 ^ [{2n - l)72„-i/i2n-i - (2n)72„/i2„] • 



(C.9) 



n>l 



Replacing /?.2n-i and /i2n by their expressions (1C.4I) and taking into account (I2.10p we obtain 
that 7i is given by the integral involving the functions 7±(t). It takes much simpler form when 
expressed in terms of the functions T±{t) defined in (12.41) 



7i 



m / 

/ dt 



■ng 



^2 _|_ ■J^2g2 



(r_(t)-r+(t)) + 



-j-2 _|_ ■J^'2.g2 



Making use of identities 



i2 + 7rV 
t 

t2 + 772^2 



we rewrite 71(5') as 



7i 



777 / 

/ dwe"^^" 



/■oo 

/ rfw e-'^s" cos {ut) , 
Jo 

/■oo 

/ du e-^'s^ sin (wt) , 
Jo 

/•oo 

/ rft cos(Mt) (r_(t) - r+(t)) 

.Jo 

dt sin(Mt) (r_(t) + r+(t)) 



(C.IO) 



(C.ll^ 



(C.12) 



Let us spit the m— integral into < m < 1 and m > 1. We observe that for u'^ < 1 the t— integrals 
in this relation are given by (12.61) . Then, we perform integration over u > 1 and find after some 
algebra (with T{t) = r+(t) + iT_{t)) 



7i 



2m' 



(1 - e-"^') 



V2^ 



m 



-■ng 



Re 



TT 



dt 



uo 



t + mg 



(C.13) 



Substituting this relation into ( 1C.7I) we arrive at (I4.36p . 



D Relation to Whittaker functions 

In this appendix we summarize properties of special functions that we encountered in our analysis. 
Integral representations 



Let us first consider the functions Vn{x) (with n = 0, 1) introduced in (I2.4ip . As follows from their 
integral representation, Vo{x) and Vi{x) are entire function on a complex x— plane. Changing the 
integration variable in (12.411) as m = 2t — 1 and -u = 1 — 2t we obtain two equivalent representations 



Vnix) = l23/2-e" /' rftr 1/4(1 _t)l/4~n, 
TT Jn 



-2tx 



(D.l) 
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which give rise to the following expressions for Vn{x) (with n = 0, 1) in terms of Whittaker 
functions of the first kind 

T(- - n) 

Vn{x) = 2"" . (2x)"/^-^M„/2-l/4.1/2-n/2(2x) , 

r(f)r(2 -n) 

= 2-" ^^5)r(2 ^n) (-2^)"^'"'^V4-n/2,i/2-n/2(-2x) . (D.2) 

In distinction with Vn{x), the Whittaker function Mn/2-i/4,i/2-n/2(2a;) is an analytical function 
of X on the complex plane with the cut along negative semi-axis. The same is true for the factor 
^2a;)"-/2-i gQ f^ii^g^f^ the product of two functions in the right-hand side of (]D.2|) is a single-valued 
analytical function in the whole complex plane. The two representations (ID. 20 are equivalent in 
virtue of the relation 

M„/2_i/4,i/2-n/2(2x) = e^'^^^^-^/^) Mi/4_„/2,i/2-n/2 (-2x) (for Imx ^ 0) , (D.3) 

where the upper and lower signs in the exponent correspond to Ima; > and Imx < 0, respec- 
tively. 

Let us know consider the functions U^{x) and U^{x). For real positive x they have an 
integral representation fl2.4ip . It is easy to see that four different integrals in fl2.4ip can be found 
as special cases of the following generic integral 

1 r°° 



2 _ 

defined for a; > 0. Changing the integration variable asu = t/a; + lwe obtain 

Uabix) = 2-+^-3/2a;-''-i/2 J dt e-* ( 1 + ^ j • (D.5) 

The integral entering this relation can be expressed in terms of Whittaker functions of second 
kind or equivalently confluent hypergeometric function of the second kind 

Uab{x) = 2^-3/2r(i - a + b)x-'~'/^ e^ Wab{2x) , (D.6) 
= ir(i - a + 6)U(| - a + 6, 1 + 26; 2x) . 

This relation can be used to analytically continue Uab{x) from x > to the whole complex 
X— plane with the cut along negative semi-axis. Matching ( ID. 41) into ( ]2.4ip we obtain the following 
relations for the functions U^{x) and U^{x) 

U+{x) = |r(|)x-^e^iy_i/4,i/2(2x), Utix) = |r(i)(2x)-i/'e^W^i/4,o(2x), 

U^{x) = |r(|)a;-ie^iyi/4,i/2(2x), U^{x) = lr{l){2x)''/^ W.y,^o{2x) . (D.7) 

The functions Vi{±x), U^{x) and Vo{±x), U^{x) satisfy the same Whittaker differential equation 
and, as a consequence, they satisfy Wronskian relations 

Vi{-x)Uo{x) - V^{-x)U{{x) = Viix)U+{x) + Vo{x)U+{x) = ^ . (D.8) 
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The same relations also follow from (IB .71) for x = ±y. In addition, 



U+{x)U^{-x) + Ut{x)U,-{-x) 



2^2; 



X 



(for Imx ^ 0) 



(D.9) 



Combining together (ID.SP and (]D.9|) we obtain the following relations between the functions 

2^/2 -^3»7r 



Vo{x) 
Viix) 



TT 

2V2 



TT 



[e^U,-(-x)+e-^U^ix)] 



(D.IO) 



where the upper and lower signs correspond to Imx > and Imx < 0, respectively. 

At first sight, the relations (1D.10I) look surprising since Vo{x) and Vi{x) are entire functions 
in the complex a;— plane, while U^{x) and U^{x) are single- valued functions in the same plane 
but with the cut along the negative semi-axis. Indeed, one can use the relations (ID.Sp and (ID.9P 
to compute the discontinuity of the these functions across the cut as 



AU^{-x) 



AU^{-x) 



IT 



±-e--'K)(Tx)^^(x) 



TT 



--e-^V,iTx)9{x) 



(D.ll) 



where AU{—x) = \im^^Q[U{—x + ie) — U{—x — ie)]/{2i) and 9{x) is a step function. Then, one 
verifies with a help of these identities that the linear combinations of t/— functions in the right- 
hand side of (ID.lOp have zero discontinuity across the cut and, therefore, they are well-defined 
in the whole complex plane. 



Asymptotic expansions 

For our purposes, we need asymptotic expansion of functions V^(x) and U^{x) at large real x. 
Let us start with the latter functions and consider a generic integral flD.6l) . 

To find asymptotic expansion of the function Uab{x) at large x, it suffices to replace the last 
factor in the integrand flD.6l) in powers of t/{2x) and integrate term by term. In this way, we 
find from flD:6|l and flDTD 



U^{x 



(2x)"^/^r(|)F(i,f|-^) = (2x)-^/^r(f) 
u,-{x) = (2x)-^/^r(|)F(-i, II - ^) = (2x)-^/^r(|) 
Utix) = (2x)-Vnr(i)F(i, il - ^) = (2x)-vnr(i) 
t/r(x) = (2x)-^/^ir(|)F(f , II - ^) = (2x)-^/^ir(|) 



32x 



+ . . . 



1 + 
1 - 
1 - 



32x 

1 

32^ 

9 
32^ 



+ . . 
+ . . . 
+ . . . 



(D.12) 



where the function F{a,b\ — ^) is defined in fl2.56p . 

Notice that the expansion coefficients in flD.12p grow factorially to higher orders but the 
series are Borel summable for x > 0. For a; < one has to distinguish the functions U^{x + ie) 
and U^{x — ie) (with e 0) which define analytical continuation of the function U^{x) to the 
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upper and lower edges of the cut, respectively. In contrast with this, the functions Vn{x) are 
well-defined on the whole real axis. Still, to make use of the relations (ID. 101) we have to specify 
the [/—functions on the cut. As an example, let us consider Vo(— tt^j) in the limit g oo and 
apply flmOj) 



2V2 



4 e 



TT 



[U+{-ng + ze)+e-'-^Uo{7cg)] 



(D.13) 



where e — > and we have chosen to define the [/—functions on the upper edge of the cut. Written 
in this form, both terms inside the square brackets are well-defined separately. Replacing 
functions in flD.lSp by their expressions flD.12p in terms of F— functions we find 



Voi-Tig) 



{2ng)-^/^ e^9 



1 3| 
4' 4I 



277g J 



(D.14) 



with given by 



A^ 



a 



r(f) 



i27rg) 



1/2 



cr 



3i7V 

e 4 



(D.15) 



Since the second term in the right-hand side of ( 1D.14[) is exponentially suppressed at large g 
we may treat it as a nonperturbative correction. Repeating the same analysis for Vi{—7Tg), we 
obtain from (fPAnll and (lDl2l) 



Vi{-ng) 



(27r^)-5/^ e^s 



2r(|) 



8ngF 



1 1 1 



3 3 1 



1 

27rg 



(D.16) 



We would like to stress that the '-|-ze' prescription in the first term in flD.14p and the phase factor 
a = e in (ID.ISP follow unambiguously from (1D.13|) . Had we defined the [/—functions on the 
lower edge of the cut, we would get the expression for Vo(— vr^f) with '— ie' prescription and the 
phase factor . The two expressions are however equivalent since discontinuity of the first term 
in flD.14p compensates the change of the phase factor in front of the second term 



1 51^ 

4' 4 I 27rg ^ ^ 



' 4' 4 I 27rg 



a 



- ^1 

45 4I 



1 

27rg 



(D.17) 



If one neglected '+ie' prescription in (]D.13|) and formally expanded the first term in (10.140 
in powers of 1/g, this would lead to non-Borel summable series. This series suffers from Borel 
ambiguity which are exponentially small for large g and produce the contribution of the same 
order as the second term in the right-hand side of (lD.14p . The relation (lD.14p suggests how to 
give a meaning to this series. Namely, one should first resum the series for negative g where it is 
Borel summable and, then, analytically continue it to the upper edge of the cut at positive g. 



E Expression for the mass gap 

In this appendix we derive the expression for the mass gap (14.41) . To this end, we replace 
r {An git) in (14.11) by its expression (14.21) and perform integration over t in the right-hand side of 
(14.11) . We recall that, in the relation (14. 2p . Vo^i{4:TTgt) are entire functions of t, while fo^i{4:7igt) 
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are meromorphic functions defined in (I4.3p . It is convenient to decompose T^Angit)/ (t + |) into 
a sum of simple poles as 



r{A7cgit) 



. . ^ V^fc(47r^t) fki4:Trgt) - fk{-ng) 
fk{-Trg) + VkiAngt) 



fc=0,l 



fc=0.1 



where the second term is regular at t = — 1/4. Substituting this relation into (14. 11) and replacing 
fk{4:7Tgt) by their expressions fl4.3p . we encounter the following integral 



Rk{4:7Tgs) = Re 



dt e-4-9*-^-/4 



t-s 



= Re 









dt e-*-*"/^- 



Vkjt) 
t — 4:7igs 



Then, the integral in (14. ip can be expressed in terms of i?— function as 



Re 



dt e 



Airgt—iTT /A 



fo{-ng)Ro{-'Kg) + fi{-Tig)Ri{-'Kg) 



n>l ^^+4 



+ Yl ^^^7^ [Ur{Arrng)R^{-A7rgn) - U, {Anng)R,{-A7rgn)] 



n>i 



(E.2) 



(E.3) 



where the last two lines correspond to the second sum in the right-hand side of (lE.ip and we 
took into account that the coefficients c±{n,g) are real. 

Let us evaluate the integral flE.2l) and choose for simplicity Ro{s). We have to distinguish 
two cases: s > and s < 0. For s > we have 



Rn(s] 



Re 
72 



/ dv e-(i~^)^ 



dt e-'^* •l/o(t) 



Re 



TT 



du 



'l + ufl\l-u)-^l^' 



u — V — le 



:e.4) 



where in the second relation we replaced Vo(t) by its integral representation fl2.4ip . Integration 
over u can be carried out with a help of identity 



1 (l+^)l/4-/c(l_^)-l/4 , 

du ^ '- ^- — '- = (5fc,o - + 1 X 



+nV4 



> 1 



In this way, we obtain from flE.4p 



/4(1±£)V4 



:e.5) 



v' < 1 



dv e-^^-^^' 



v + 1 
v-1 



1/4 



V2 



-2e-2^ mis) 



(E.6) 



with the function Uq{s) defined in fl2.4ip . In the similar manner, for s < we get 



Rois) V2 



- + 2U,-{-s) 
s 



(E.7) 
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together with 

Ri{s) = 2V2 [e{-s)U^{-s) + e{s) e-2- U^{s)] . 
Then, we substitute the relations (lE.6p . (IE. 70 and (lE.Sp into (1E.3P and find 

^r(47r(72t)" 



(E.8) 



Re 



2y/2M-ng) 



2ng 



+ 2V2 f,{-ng)U{{ng) + ^ [/o(-7r(7) + 1] 

ng 



:e.9) 



where the last term in the right-hand side corresponds to the last two lines in (1E.3P (see 
Eq. fl2.4Up ). Substitution of (1E.9P into (14. ip yields the expression for the mass scale fl4.4p . 
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